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CHAPTER  I 


INTRODUCTION 


Problems  related  to  liquid  mixing  in  various  systems  have  attracted  the 
interest  of  many  researchers.   While  the  importance  of  the  mixing  mechanism 
has  been  recognized,  knowledge  of  the  mechanism  has  not  been  extensively 
developed.   There  has  been  great  impetus  for  further  study  of  liquid  mixing, 
because  a  thorough  understanding  of  the  mixing  mechanism  within  a  system  of 
mass  transfer  or  chemical  reaction  is  helpful  in  successful  design,  operation* 
and  control.   Liquid  mixing  on  distillation  trays  is  so  important  that  neg- 
lecting it  often  brings  about  an  appreciable  error  in  the  prediction  of 
Murphree  efficiency  and  thus  leads  to  poor  design  of  the  distillation  column. 

This  work,  a  continuation  of  previous  work  (1),  presents  results  of  a 
study  of  liquid  mixing  on  a  non- ideal  sieve  tray.   The  term  "non- ideal  sieve 
tray"  is  used  here  to  distinguish  it  from  a  commercial  sieve  tray.   This  tray 
has  fewer  holes  or  a  smaller  percentage  of  free  area  than  commercial  ones. 
The  use  of  such  a  tray  may  enable  one  to  develop  a  clearer  physical  picture 
of  fundamental  mechanism  of  flow  patterns. 

Since  Danckwerts  (2)  introduced  the  concept  of  age  distribution  analysis 
to  the  study  of  flow  systems,  much  work  has  been  based  on  this  concept. 
Danckwerts  originated  the  use  of  two  probability  functions,  the  exit  age  dis- 
tribution and  the  internal  age  distribution  function.   The  present  study  of 
liquid  mixing  on  a  distillation  tray  uses  these  two  functions  which  are  dis- 
cussed in  detail  in  Chapter  3.   Noar  and  Shinnar  (3)  employed  the  intensity 
function,  which  represents  another  approach  to  the  study  of  fluid  mixing. 

Due  to  the  complexity  of  the  distillation  tray  system,  an  analytic 


solution  of  the  governing  equation  for  liquid  mixing  in  the  system  is  very 
difficult,  if  not  impossible,  to  obtain*   The  approximate  solution  is  based 
on  the  fact  that  the  extent  of  actual  mixing  lies  between  two  extreme  modes. 
These  two  extreme  modes  are  represented  by  the  two  models:  the  plug  flow 
model  and  the  backmix  (completely  stirred)  model.   The  combined  model  which 
describes  the  actual  mixing  is  made  by  series,  parallel,  or  series-parallel 
combination  of  these  two  extreme  models  and  is  usually  modified,  for  instance, 
by  a  bypassing  stream  or  a  dead  region.   The  proposed  models  are  developed 
in  the  above  way  to  fit  the  experimental  data. 

The  overall  intent  of  this  work  was  to  lay  the  foundation  for  a  com- 
prehensive investigation  of  liquid  mixing  on  the  non-ideal  sieve  tray.   The 
specific  purposes  of  the  work  carried  out  here  were  the  following: 

(1)  To  develop  an  experimental  procedure  for  the  pulse  testing  tech- 
nique by  using  the  conductivity  method. 

(2)  To  develop  flow  patterns  in  order  to  obtain  a  qualitative  picture 
of  their  effects  on  liquid  mixing  by  means  of  a  photographic  technique  and 
visual  observation. 

(3)  To  detect  the  possible  existence  of  the  stratification  effect  in 
the  system  and  the  extent  to  which  it  may  disappear. 

(4)  To  suggest  a  technique  or  new  method  of  monitoring  a  tracer  in 
case  the  stratification  does  exist. 

The  pulse  testing  technique  measures  mainly  the  concentration- time 
curves  at  the  inlet  and  at  the  outlet  of  the  system.  The  concentrations 
across  the  inlet  and  the  outlet  may  be  non-uniform  due  to  the  stratified 
streams  which  occur  in  the  system.  The  nonuniformity  of  concentrations 
across  any  cross-sections  may  affect  the  degree  of  mixing  and  model 


parameters.   For  this  reason,  the  effect  of  those  stratified  streams  or  the 
effect  of  stratification  (bulk  segregation)  on  the  degree  of  liquid  mixing 
needs  to  be  investigated.   The  detection  of  the  stratification  effect  in  the 
system  was  done  qualitatively  by  photography  and  by  visual  observation.   A 
quantitative  measurement  of  this  effect  has  been  carried  out  experimentally. 
In  general,  the  major  accomplishment  of  this  research  was  the  establish' 
ment  of  the  equipment  and  instrumentation,  followed  by  a  primary  detection 
of  stratification  effect;  a  minor  contribution  was  a  statistically  designed 
plan  to  obtain  experimentally  the  model  parameters* 


CHAPTER  II 

STRATIFICATION  EFFECT 

As  has  been  pointed  out  by  Bell  and  Babb  (4),  the  effect  of  flow  strati- 
fication is  important  in  monitoring  tracer  experiments*   There  are  two  ways 
of  monitoring  tracer  experiments.  One  of  these  is  to  place  a  detector  either 
at  the  wall  or  in  the  inner  part  of  the  vessel  in  such  a  way  that  it  monitors 
the  average  or  mixing  cup  concentration  of  the  vessel  contents  in  some  plane 
through  the  vessel.   The  other  is  to  measure  the  mixing  cup  concentration  of 
the  effluent  leaving  a  plane  through  the  vessel.   The  first  method  was  called 
"through  the  wall"  and  the  second  "in  the  effluent".   Let  one  consider  that 
the  flow  is  made  up  of  many  stratified  elements.  Assuming  that  a  pulse  of 
tracer  has  been  injected,  each  fluid  element  contains  a  certain  quantity  of 
that  initial  injection.  The  structure  of  the  flow  may  appear  as  shown  in 
Fig.  1  which  represents  a  plane  cut  through  the  weir  at  the  detection  posi- 
tions. The  stratified  elements  are  shown  in  Fig.  1,  the  length  of  each  of 
the  stratified  elements  being  arbitrary.  The  velocities  of  the  elements,  v., 
the  concentrations,  c,,  and  the  areas,  a,,  may  all  differ.  The  volumetric 
flow  rate  Q  can  be  represented  by 

E.  /     v  da   =  Q  (2-1) 

a.(t)   l 

i 

If  the  concentration  is  monitored  through  the  wall  as  a  function  of 
time,  the  contributions  of  each  element  to  the  monitored  mean  concentration 
are  the  cross -sectional  area  and  concentration  of  the  element.  By  denoting 
this  mean  concentration  by  C.(t),  one  can  write 


X    Detection   position 


wall 


Fig.  I.  Schematic  representation  of 
stratified   flow   of  tracer 
elements   in  a  plane  , 


Z'/     C.(t)da 

i-.t<t) 


C.(c)  =  =-2/  C,(t)da  (2-2) 


E  J     da        A  i  a 
i  a^t) 


i 


where  A  denotes  the  cross-sectional  area* 

If  the  concentration  is  monitored  in  the  effluent  also  as  a  time  function, 
the  contributions  of  each  element  to  the  mean  concentration  are  the  area,  con- 
centration, and  the  velocity  of  the  element.  By  denoting  this  mean  concen- 
tration by  C„(t),  one  has 

Z  J  c,<t)v.  da 

C2(t)  =  —  ="£/     C.(t)v.da    (2-3) 

S  /    v.  da       Q  i  a,  (t)  *    * 
iai(t)  *  l 

It  can  be  seen  that  the  two  monitoring  techniques  must  be  identical  in 
the  case  that  there  is  no  concentration  gradient  or  no  velocity  gradient 
across  the  cross-section,  and  thus  C. (t)  =  C«(t).  Actual  situations  in  liquid 
mixing  do  not  correspond  to  these  two  extreme  conditions  of  no  concentration 
gradient  and  no  velocity  gradient*   Therefore,  the  effect  of  stratification 
or  bulk-segregation  on  the  degree  of  liquid  mixing  needs  to  be  investigated. 

Let  us  consider  the  detection  of  the  stratification  effect  on  a  specific 
system  in  which  pulse  testing  is  carried  out  with  the  use  of  two  probes,  one 
located  upstream  and  the  other  downstream.  To  employ  the  pulse  testing 
technique  to  measure  the  stratification  effect  only  two  probes  have  been  used 
simultaneously  in  this  work  because: 

(1)  The  stratification  may  occur  either  upstream  or  downstream, 

(2)  A  small  number  of  probes  minimizes  disturbance  of  the  flow  pattern 
by  the  probes  being  used. 


CHAPTER  III 


AGE  DISTRIBUTION 


Knowledge  of  the  complete  flow  pattern  of  fluid  passing  through  a  vessel 
can  be  obtained  in  different  ways.  A  straightforward  way  would  be  to  tag  and 
follow  each  individual  molecule  as  it  passes  through  the  vessel.   In  this 
approach  a  complete  velocity  distribution  of  the  fluid  within  a  vessel  could 
be  obtained.   However,  the  attendant  complexities  make  it  impossible  to  use 
this  approach. 

An  alternate  approach  is  based  on  the  concept  of  age  distribution  ana- 
lysis (2).   This  concept  can  provide  information  concerning  the  fraction  of 
fluid  that  resides  for  a  certain  time  in  the  vessel.   The  details  of  the 
exact  path  history  or  where  the  particle  has  been  during  its  stay  are  not 
considered.   Therefore,  information  about  the  point  to  point  changes  of  the 
variables  is  not  obtainable  from  this  approach.   However,  the  concept  of  age 
distribution  can  provide  us  with  the  exit  age  distribution  function.  A 
mathematical  model  based  on  this  function  and  related  to  the  physical  situ- 
ation is  then  set  up  to  fit  experimental  data. 

Due  to  the  linearity  of  the  stimulus-response  of  the  mixing  behavior, 
an  injected  tracer  (assuming  no  absorption  at  the  walls  and  no  disappearance 
by  reaction)  will  pass  through  the  vessel  with  the  rest  of  the  fluid.   The 
age  distribution  of  the  tracer,  therefore,  represents  the  distribution  of 
the  fluid  elements  of  the  system.  Due  to  this  property,  the  age  distribution 
can  be  used  to  understand  the  flow  pattern. 

It  is  necessary  to  define  the  mean  residence  time  and  the  reduced  time 
before  proceeding  to  the  discussion  of  the  age  distribution  (3)« 
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Definition  of  the  mean  residence  time 

The  mean  residence  time,  t,  is  defined  as 


t  ■  -  —  -       Volume  of  the  vessel  available  for  flow  ,~    . k 

v     Volumetric  flow  rate  of  fluid  through  the  vessel 


The  reduced  time,  0,  is  defined  as 


0=4-  (3-2) 

t 


The  discussion  following  is  based  on  the  work  of  Levenspiel  <8)» 

Internal  age  distribution  (2) 

A  vessel  contains,  in  general,  fluid  elements  of  varying  ages,  age  being 
the  time  a  fluid  element  has  spent  in  the  vessel.   Let  the  symbol,  .1(8),  be 
the  measure  of  the  distribution  of  the  ages  of  the  fluid  in  the  vessel  defined 
in  such  a  manner  that  J_(0)  d0  is  the  fraction  of  fluid  elements  in  the  vessel 
between  the  ages  0  and  0  +  d0.  A  typical  plot  of  the  I(Q)   against  6  is  given 
in  Fig.  2.   Since  the  sum  of  all  the  fractions  of  the  fluid  elements  in  the 
vessel  is  unity,  this  sum  must  also  be  the  total  area  under  the  ,1(9)  vs.  0 
curve.   Thus 

oo 

/  1(9)  d9   =   1  (3-3) 

0 

The  fraction  of  the  fluid  elements  with  the  ages  younger  than  a  certain  value, 
say  0  ,  is  shown  in  Fig.  2  as  the  shaded  area  and  is  given  by 

e. 

;  l  ice) de  (3-4) 

o 


i(e) 


Dimensionless    Time,  0 

Fig. 2.  A    typical  internal  age    distribution 
curve. 


E(e) 


D!.7icr:c;or;]QC3    TirnQ,  © 


iQ.3.  a   typical    e^jt    age   distribution 
curve . 
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The  fraction  of  the  fluid  elements  with  ages  older  than  0 .  is 


/  1(9)  d6   =  1  -  /  l   1(9)  d0  (3-5) 

e  o 

The  internal  age  distribution  I(t)  based  on  actual  time  rather  than  reduced 

time  is  related  to  1,(6)  as 

i 

1(9)  =  tl(t)       with     /  I(t)dt   a  i         (3-6) 

0 


Exit  age  distribution  (2) 

In  a  manner  similar  to  the  internal  age  distribution  function,  let  E(0) 
denote  the  distribution  of  age  of  all  fluid  elements  leaving  the  vessel* 
E(9)  is  defined  in  such  a  manner  that  E(9)d0  is  the  fraction  of  fluid  ele- 
ments in  the  exit  stream  between  the  ages  of  0  and  0  +  d0.   It  follows  that 


oo 

/  E(0)d0   =   1  (3-7) 

0 


A  typical  E(0)  vs.  0  curve  is  given  in  Fig,  3,  The  fraction  of  fluid 
elements  in  the  exit  stream  with  ages  younger  than  age  02  is 

e2 

/ZE(0)d0  (3-8) 

0 

The  fraction  of  fluid  elements  with  ages  older  than  92  is  then 


/  E(8)<ie   =  1  -  /   E(e)d6  (3-9) 

62  0 
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E(9)  Is  variously  referred  to  as  the  exit  age  distribution  function  or  the 
exit  residence  time  distribution  or  simply  the  residence  time  distribution 
function  (r.  t.  d.  f.).  If  actual  time  t  is  used  instead  of  d linens lonless 
or  reduced  time  6, 

00 

E(0)   =  tE(t)       with      /  E(t)dt   =   1 

0 

Intensity  function  (3) 

Another  method  for  representing  residence  time  variability  In  continuous 
flow  systems  is  based  on  the  use  of  the  intensity  function  (3).   The  advantage 
of  this  method  is  that  it  allows  a  physical  insight  into  the  mixing  processes 
within  the  system  and  enhances  the  interpretability  of  experimental  curves. 

The  intensity  function  ^(9)  is  defined  in  such  a  manner  that  ^(9)d0 
represents  the  fraction  of  fluid  elements  in  a  vessel  of  age  9  that  will 
leave  at  a  time  between  0  and  9  +  d9. 

The  physical  meaning  of  the  three  probability  functions 

The  amount  of  material  leaving  in  time  less  than  0  is 

9 
v/  E(9)d9  (3-10) 

0 

where  v  =  volumetric  flow  rate,   Thus,  the  amount  of  fluid  not  leaving  the 
vessel  (i.e.  remaining  in  vessel)  is 

9 
vl(9)   »  v-  v/  E(0)d9  (3-11) 

0 

and  thus 


12 


e 
1(e)    =    1  -  /   E(e)de  (3-i2) 

o 


Consequently,  one  can  write 


><e>   =  -2^1  (3-13) 

-       de 


The  amount  of  material  leaving  the  vessel  between  0  and  9  +  d9  can  be 
written  as  a  product  of  the  amount  not  leaving  before  9  and  the  fraction  of 
fluid  of  age  6  that  will  leave  between  9  and  9  +  d9.   In  other  words, 

vE(9)d9  =  [vl(9)]f\(9)  d9] 
or  (3-14) 

\(9)   = 


E(9)        t   dl(9)      din  1(9) 
K9)  "~   '  1(e)  d9  d9 


The  F  curve  (2) 

With  no  tracer  initially  present  in  the  entering  stream,  let  a  step  of 
tracer  signal  of  concentration  C  be  introduced  into  the  fluid  entering  the 
system  in  such  a  manner  that  the  volumetric  flow  rate  to  the  vessel  remains 
constant.  Then  the  concentration-time  curve  for  the  tracer  in  the  exit  fluid 
stream,  measuring  tracer  concentration  in  terms  of  inlet  tracer  concentration 
and  measuring  time  in  reduced  units,  is  called  the  F-curvee  As  shown  in 
Fig.  4  the  F-curve  rises  from  0  to  U 

The  C  curve  (2) 

The  curve  that  describes  the  concentration- time  function  of  tracer  in 
the  exit  stream  of  a  vessel  in  response  to  an  idealized  delta  function  of 
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Step  function 


II 


tracer  input  signal 

/      Tracer   output 
/,        signal   or  Fcurve 

Dimensionless    timef0 

Fig.  4.   Typical   downstream  signal,    the    £•    . 
curve,     in     response  to  an  upstream 
step    input   signal. 


Delta    function  or  pulse 
tracer    input    signal 


Tracer  output  signal 
or    C    curve 


Dimensionless     time  ,  0 


Fig.  5.    Typicd     downstream    signal,     the 

£  curve,   ['in        response    to  an   upstream 
delta -function    tracer    input     signal. 
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unit  impulse  input  is  called  the  C-curve  if,  as  in  the  case  of  the  F  curve, 
dimensionless  coordinates  are  chosen.  Concentrations  are  measured  in  terms 
of  the  initial  concentration  of  injected  tracer,  C  ,  as  if  it  were  evenly 
distributed  throughout  the  vessel,  and  time  is  measured  in  reduced  units. 
The  area  under  the  C-curve  always  be  unity,  that  is, 


00  OO  00  00 

/  C(0)de  =  /  -7-dO  =   1  ,   where   C  =  £  /  Q(t)dt  =  4-/  C(t)dt 
0  -         0  Co  °    V0  tO 


Relation  between  the  F,  C,  J[,  and  E  curves  in  closed  vessels  (8) 

By  material  balance  the  following  relationships  can  be  obtained* 

F(9)  +  1(9)  =  1 


c<„  „  ,(,-,  .  S<!>  =  .JS2>  (3-i5) 

~        d0         d0 


e  e 

F(9)  =  1-1(9)  =  /  E(9)d0   =  /  C(0)d0 

0  0 
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CHAPTER  IV 
SYSTEM  AND  TRACER  TECHNIQUE 

Typo  of  flow 

There  are  two  extreme  modes  of  mixing,  micromixing  and  macromixing.   The 
former  is  mixing  on  the  microscopic  scale  (mixing  of  individual  molecules), 
while  the  latter  is  mixing  on  the  macroscopic  scale  (mixing  of  clumps,  groups, 
or  aggregates  of  molecules).   An  actual  mixing  condition  of  a  fluid  usually 
is  found  to  be  between  the  two  extreme  cases  and  is  called  a  partially  segre- 
gated fluid.  A  microfluid  exhibits  no  segregation  and  a  macrofluid,  hence, 
exhibites  complete  segregation.   The  degree  of  segregation  is  of  importance 
in  analyzing  the  actual  conversion  of  the  reactions  in  nonlinear  rate  pro- 
cesses (5).   However,  while  most  of  the  mass  transfer  processes  are  linear 
processes,  mixing  on  a  distillation  plate  can  be  treated  from  either  point 
of  view  (9). 

The  pulse  testing  technique  used  in  detecting  the  residence  time  distri- 
bution has  been  extended  to  a  two  phase  and  even  to  a  three  phase  system 
(two  phase  flow  in  fluidized  beds)  (10).   The  treatment  of  residence  time 
distribution  in  a  two  phase  system  (for  instance,  a  distillation  plate)  can 
be  simplified  by  independent  treatment  of  each  individual  phase  because  the 
mixing  process  is  a  linear  process.   The  first  phase  (the  gas  phase)  offers 
an  additional  source  of  kinetic  energy  as  well  as  an  additional  impediment 
to  the  flow  of  the  second  phase  (the  liquid  phase)  (11)- 

Open  and  closed  systems 

An  open  system  is  essentially  a  system  with  no  discontinuity  in  the 
type  of  flow  at  the  point  of  tracer  injection  or  at  the  point  of  tracer 
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measurement,  i.e.  diffusion  and  dispersion  have  to  be  taken  into  consideration 
in  addition  to  bulk  flow  at  the  points  of  injection  and  measurement.  On  the 
other  hand,  a  closed  system  is  defined  as  one  in  which  a  fluid  enters  and 
leaves  by  bulk  flow  alone.  When  either  the  entrance  or  exit  alone  satisfies 
the  closed  system  requirements,  the  system  is  said  to  be  closed-open  or  open- 
closed. 

A  distillation  system  can  be  considered  to  be  a  closed  system.  This  is 
based  on  the  assumption  that  diffusion  and  dispersion  are  negligibly  small 
at  the  downcomers  as  compared  with  that  on  the  active  zone  of  the  tray. 
Therefore,  the  transport  of  the  tracer  at  both  entrance  and  exit  can  be  con- 
sidered to  be  due  to  the  bulk  flow  only. 

Tracer  technique 

The  use  of  tracer  experiments  for  characterizing  flow  patterns  in  process 
equipment  has  become  a  subject  of  much  interest.   The  system  is  perturbed  with 
the  different  kinds  of  tracer  (dye  stuffs,  electrolytes,  and  radioactive 
isotopes),  and  the  concentration- time  curve  at  the  exit  point  or  at  both 
entrance  and  exit  points  is  analyzed,,  There  are  several  methods  extensively 
used  for  the  system  excitation:  step  input  of  tracer  (or  heaviside  step  func- 
tion input),  delta  function  input  of  tracer,  sinusoidal  input  of  tracer,  and 
arbitrary  input  of  tracer. 

The  first  two  methods  of  tracer  input  for  the  system  excitation  are 
directly  related  to  the  responses  which  give  rise  to  the  F-curve  and  the 
C-curve  shown  in  Figs.  4,  5.   That  is  to  say,  the  F-curve  and  the  C-curve 
can  be  obtained  by  measuring  only  the  concentration- time  curves  of  the  output 
responses.   However,  the  techniques  using  the  definite  shapes  of  tracer  input 
signals  have  their  difficulties  and  inaccuracies.  While  the  sinusoidal  input 
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of  tracer  is  usually  used  for  the  frequency  response  analysis,  there  are  also 
several  difficulties  and  disadvantages  associated  with  the  method  (1,  12,  13, 
14). 

Pulse  testing  technique 

In  the  pulse  testing  method,  the  tracer  is  injected  in  a  very  arbitrarily 
manner  at  a  convenient  point  upstream  from  the  system  of  interest.   Since  the 
input  is  arbitrary,  an  arbitrary  response  signal  is  obtained,  and  thus  the 
measurement  of  the  concent rat ion- time  curve  is  no  longer  limited  to  the  outlet 
of  the  system  only.   In  many  situations  the  measuring  points  are  located  at 
the  inlet  and  at  the  outlet,  and  the  measured  curves  are  denoted  as  X(t)  and 
Y(t)  respectively.   Their  shapes  are  apparently  arbitrary,  but  their  magni- 
tudes are  greater  than  zero  during  the  measuring  period.   Typical  input  and 
output  concentration- time  curves  are  shown  in  Fig.  6. 

The  dimensionless  variables  corresponding  to  X(t)  and  Y(t)  can  be  defined 
as 

x<e)  =  ^2 

A 
x 


where 


Y(9)  = 


t  =  I 
v 


e   =  4 


(4-1) 

tY(t) 


and 


A   =  /  X(t)dt 


0 


A   =  /  Y(t)dt 
y     0 
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CHAPTER  V 
THE  MOMENTS  METHOD  OF  ANALYSIS 

The  interesting  properties  of  the  moments  method  of  analysis  have  been 
pointed  out  by  several  authors  (10,  15).   This  method  is  particularly  well 
adapted  to  the  description  of  the  quality  of  mixing  in  a  continuous  mixing 
vessel.   It  is  possible  to  describe  numerically  the  process  dynamics  without 
the  necessity  of  a  specific  model.   On  the  other  hand,  if  a  system  model  is 
presumed,  the  system  moments  can  be  used  to  determine  coefficients  without 
the  use  of  curve  fitting  techniques.   The  characteristics  of  the  residence 
time  distribution  function  can  be  represented  by  a  set  of  moments  and  con- 
sequently by  a  set  of  numbers. 

The  first  moment,  or  mean,  of  the  residence  time  distribution  function, 
E(9),  is  given  by 

00 

w    =  Mj  =  /  9E(9)d9  (5-1) 

The  second  moment  about  the  mean,  commonly  called  the  variance,  measures 
the  spread  of  the  distribution  about  the  mean  and  is  equivalent  to  the  square 
of  the  radius  of  gyration  of  the  distribution.   It  is  defined  for  a  continuous 
distribution  as 

G2Q     =  M2   =   /   (9  -yUe)2  E(0)d9  (5-2) 

The  third  moment  about  the  mean,  commonly  called  the  skewness,  is 
obtained  as 

XmQ     =  M3  =  /  (9  y/e)3E(9)d9  (5-3) 
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The  exit  age  distribution  function,  E(9),  is  a  mathematical  representa- 
tion of  a  model  derived  from  a  physical  system.   In  the  model  there  exist 
some  parameters  which  we  have  to  be  determined  experimentally.   In  order  to 
determine  the  model  parameters  and  fit  the  distribution  curve  to  the  experi- 
mentally collected  input  and  output  data,  X(9)  and  Y(9),  the  transfer  function, 
H(s),  of  a  mixing  system  has  to  be  introduced.   It  is  defined  as 

H(s)  =  =^4  (5-4) 

X(s) 


where 


-  -s9 

Y(s)  =  Laplace  transform  of  Y(9)  =  /  Y(9)e  '  d9 

0 


-  °°      -s9 

X(s)  =  Laplace  transform  of  X(9)  =  /  X(9)e    d9 

0 


If  X(9)  is  the  delta  function,  then 

X(s)  =  1  (5-5) 

and 

Y(s)  =  C(s)  =  E(s)  (5-6) 

Thus  it  is  evident  that 

E(s)  =  H(s)  (5-7) 

that  is,  the  Laplace  transform  of  the  residence  time  distribution  function 
is  equal  to  the  transfer  function  of  the  system.   From  this  identity  the 
following  relationships  will  be  developed s 

(1)  The  relationship  between  the  residence  time  distribution  function 
and  the  moments* 
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(2)  The  relationship  between  the  transfer  function  and  the  moments. 

(3)  The  relationship  between  the  moments  from  case  1  and  from  case  2. 
One  defines  the  Laplace  transform  of  the  residence  time  distribution 

function  as 

E(s)  =  /  E(9)e'S  d0  (5-8) 

0 

Upon  successive  differentiation  of  Equation  (5-8),  one  has 


dE(s)        °°  -     _  -s9 

-= =  -  /  [0E(9)]e   d9  (5-9a) 

ds         A 


d2E(s)       °°  _  2     -,  -s9 

"2         »  +/  [0  E(9)]e    d9  (5-9b) 


and 


ds         0 


.3-, 


i^i  =  -  /  [93E(9)le-S9d0   '  (5-9c) 

dsJ        0 


Taking  the  limit  on  both  sides  of  Equations  (5-9a),  (5-9b)  and  (5- 9c),  one 
obtains 

dE(s) 

lim~ =  -  /  [9E(9)]d9  (5-lOa) 

s-0  ds        o 

2- 
d  E(s)       °°  _  2 

lim  — = =  +  /  [0  E(9)ld9  (5-lOb) 

s-0   ds        0 


and 


d  E(s)       °°   3 

lim  _= =  -  /  [9  E(9)]d9  (5-lOc) 

s-^0   ds        o 
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Substituting  Equation  (5- 10a)  into  Equation  (5-1)  yields 

dE(s) 

U   =  -  lim  -= (5-lla) 

r   ©      s+0     ds 

Substituting  Equations  (5- 10a)  and  (5-l0b)  into  Equation  (5-2)  and  making 


2 


0  =  /  (6  -yUe)  E(0)d0 


00     o  oo  0     °° 

=     /     6    E(0)d0  -  2/LJ     9E(e)d0  +  /,,.  /     E(0)d9 
0  /00  /eo    ~ 

,,      rd  I<s>        o,dI<s\2        ,d£(s),2 

*     lim  I — =^—    -  2(-= )     +  (-= )   •  l  | 

s-»0         ds2  ds  ds 

» 

2-  - 

nd  E(s)  dE(s)  2-, 

=     lim  [     -\'  -    <-=i-ir]  .  (5-llb) 

s+0         ds  ds 


Similarly,  making  use  of  Equations    (5-3),    (5-1),    (5-2),    (5-l0a),   (5- 10b),  and 
(5- 10c)  yields 

3  ■         rd3E(s)         odE(s)d2E(s)        0  d|(s)   3-.  /B      ,    x 

7-        =     -  lim  [     "  -  3-=- =_  +  2(-= )].  (5-llc) 

Ue  s*0  L     ds3  ds         dg2  ds 

Now,   let  us  establish  the  relationship  between  the  transfer  function  and 
the  moments. 

Differentiation  of  Equation  (5-4)  yields 


dH(s)     =     X(s)Y'(s)   -  Y(s)X«(s) 
ds  X(s)2 
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1  Y  (      \ 

Multiplying  both  sides  of  Equation  (5-12)  by and  _     respectively 

H(s)       Y(s) 

yields 


(5-13) 


H' (s)     Y'  (s)   X*(s) 
H(s)      Y(s)    X(s) 

Upon  successive  differentiation  of  Equation  (5-13),  one  obtains 


H»(s)   H«(s)      _Y"(s)    Y'(s)  9_    rX»(s)    X»  (s)  2-,       ,e  /x 

■= — 5-  =  r (*= )   -  L- C"= )        (5-14) 

H(s)    H(s)2       Y(s)     Y(s)   J    L  X(s)     X(s) 


and 


H»»(s)   3H"(s)H«(s)   2H»(s) 
H(s)        H(sr       H(s)3 


pY">(s)  _  3  Y"(s)Y'(s)  +  2(Y«(s))3-, 
Y(s)        Y(s)2        Y(s) 


^W  _  3  X»»(s)X»(s)  +  2(XMs))3  tf 

LX(s)        X(s)2        X(s) 


In  order  to  simplify  the  above  equations,  one  introduces  the  relations 


00      -s9       °° 
limX(s)  =  lim  /  X(0)e    d9  ■  /  X(9)d9  ■  1        (5-16) 

s-»0         s-»0  0  0 


—  °°      -s9       °° 

lim  Y(s)  =  lim  /  Y(9)e  '  d9  ■  /  Y(9)d9  =  1        (5-17) 
s-0         s-0  0  0 


Hence 


lim  H(s)  =  1  (5-18) 

s-K) 
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By  taking  the  limits  of  Equations  (5-13),  (5-14),  and  (5-15)  and  making 
use  of  Equations  (5-16),  (5-17),  and  (5-18),  one  has 

lta  ^1  m  llm  pw  .  fiJsij  (5.,9) 

s-*0     ds  s-»0       ds  ds     " 

lta  [<%.  .  <Sl<£i)2]    „    lta  [^£>  .  (fiJsl)'] 

s*Q    dsz      ds         s-*0   dsz      ds 

'.lim^H£l.  (f^)2]         (5.20) 
s-0  L  ds2      ds 

ft   rd3H(s)   ,  dH(s)  d2H(s)   0,dH(g).3-, 

1  im  [ 7-  -  3 r-  +  2  (—- )  J 

s*0    dsJ      ds     dsz      ds 

s-»0    ds"3      ds     dsz      ds 

.  lta  £*<•>  .  3  «W  d^s)  +  ,^,3-j       (5.21) 

s-0  L  ds3      ds     ds2      ds 


X(s)  and  Y(s)  are  the  distribution  functions  which,  in  the  mathematical 
sense,  have  the  same  properties  as  the  residence  time  distribution  function, 
E(s).   If  one  replaces  the  E(s)  in  the  right  hand  sides  of  Equations  (5- 11a), 
(5- lib),  and  (5- lie)  by  Y(s)  or  X(s),  the  following  equations  will  be  obtained: 


u         =  -  lim  dY(s) 
/  0Y      s-*0  ds 


U  m  um  dX(s) 

•  Q*  s^o  ds 


®Y     s-K)  L  ds2     ds 
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V  ■  "3  [*% "  ^^ 

yx  s-*0         ds  ds 

(5-22) 

zl  m     .Um[^.3^^^<f>+2(^il)3l 
0Y  s-»0  ds3  ds  ds2  ds 

3  ,,      rd3X(s)        ,  dX(s)  d2X(s)        0/dX(s)   3 

ex  s-0  ds3  ds  ds2  ds 


Substituting  Equation   (5-22)    into  Equations    (5-19),    (5-20),   and    (5-21) 
yields 


i«      fd   H(s)        /dH(s),2  2  2  2 

lin  [= 5-  -   (— — )   J     =     oov  -  aov     -  A  oQW  (5-24) 


2- 
s-0  "     ds2 

and 


ii      r^3H(s)        o  dH(s)  djUs)        9  ,dH(s)   3  3  3  3  ... 

"  s'S  ^7?"  '  3  ~s 7F  +    (~}  ]   =  rev  ■  *ex   ■-«  evx        <5-25) 


Recall  that  E(s)   =  H(s).      From  Equation  (5-lla),    (5-llb),   and    (5-llc), 
it  can  then  be  established  that 

dE(s) 

-  US  "£-  = Vevx  ■  A  (5-26) 

2-  - 

d  E(s)  dE(s)   2_  2  2 

nn  [~^r  ■  c-f — >  3    =  A  <W    =    ^  <5-27> 

s-0         ds2  ds  J  9YX  0 

and 

Bl      rd  E(s)  dE(s)   d2E(s)  dE(s)   3.  3  3  ooX 

-   lim  [     ~         -  3  — =—  +  2(— )   ]     =Atftw     =  Tft  <5'28> 

s-0  L     ds3  ds  ds2  ds  0YX  ^e 
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From  the  pulse  testing  data,  Y(9)  and  X(9),  one  can  calculate  the 

2  3 

variance,  a  <JfiYX  and  the  skewness,  A  tflw'   For  a  two-parameter  model,  one 

can  make  use  of  Equations  (5-27)  and  (5-28)  to  calculate  the  model  para- 
meters. One  can  see  that,  in  general,  n+1  moments  are  required  for  a 

n- parameter  model. 

Discussion  of  the  moments  method  of  analysis 

In  general,  the  nth  moment  of  the  transfer  function  can  be  obtained  by 
taking  the  limit  of  the  nth  derivative  of  the  Laplace  transform  and  letting 
the  transform  variable  tend  toward  zero.  This  yields 


— n  °° 

lim  Y  (s)  =  (-l)n  /  tnf(t)dt 
s-0  0 


It  has  been  observed  that  the  tails  of  Y(0)  and  X(0)  curves  may  be 
unreliable  due  to  the  difficulties  of  accurate  measurements  of  tracer  con- 
centration at  low  concentrations.  Therefore,  it  becomes  obvious  that  the 
moments  method  has  its  disadvantages  (16,  17),  especially  when  higher  moments 
are  required.  Small  uncertainties  in  the  values  on  the  tail  of  the  input  and 
output  concentration  time  curves  introduce  appreciable  truncation  errors  in 
the  analysis  of  data. 

Although  the  moments  method  of  analysis  has  significant  drawbacks,  it 
is  still  of  value  due  to  a  saving  in  computer  time  in  obtaining  the  first 
and  the  second  moments  which  can  be  used  to  describe  the  extent  of  mixing* 
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CHAPTER  VI 
MODELS  OF  DISTILLATION  TRAYS 

The  models  used  in  fitting  a  mixing  system  are  of  two  kinds:  the  purely 
mathematical  model  and  the  physical  model. 

As  mentioned  by  Adler  (18),  the  use  of  a  model  has  the  following  advan- 
tages over  the  use  of  experimental  data  without  a  model. 

(1)  Scale  up  and  extrapolation  can  be  more  accurately  carried  out.   The 
model  helps  predict  the  effects  of  unavoidable  differences  in  fluid  mixing, 
which  occur  when  moving  from  laboratory  to  pilot  plant  to  a  commercial  unit* 

(2)  Experimental  data  can  be  more  fully  utilized. 

(3)  Dynamic  as  well  as  static  analysis  and  prediction  are  possible, 
permitting  early  opportunities  for  optimization  and  control  studies. 

(4)  The  amount  of  pilot  plant  time  and  data  needed  is  often  reduced. 
The  pool  model  of  a  distillation  tray,  originally  proposed  by  Kirchbaum 

(19)  and  later  by  Nord  (20),  Gatreaux  and  0,Connell  (21),  Mickley,  et  al.  (22), 
and  Leonard  (11),  postulates  that  the  plate  may  be  pictured  as  a  series  of 
completely  mixed  pools  along  the  length  of  the  tray.   While  it  assumes  com- 
plete mixing  within  each  pool,  there  is  no  mixing  between  pools.  A  plate 
with  a  single  pool  corresponds  to  a  perfectly  mixed  plate  and  one  with  an 
infinite  number  of  pools  to  a  plug  flow  plate. 

The  recycle  model  of  Oliver  and  Watson  (23)  and  Warzel  (24)  assumes  that 
the  liquid  mixing  is  affected  by  a  recycle  of  liquid  which  moves  from  the 
tray  exit  to  tray  inlet.   The  parameter  which  characterizes  mixing  is  the 
concentration  jump  at  the  inlet  weir. 

The  splashing  model  was  proposed  by  Johnson  and  Maranjozis  (25).   In  the 
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model,  splashing  of  the  liquid  is  considered  to  be  the  major  mixing  mechanism* 

The  axial  dispersed  plug  flow  model  has  been  applied  to  distillation 
plates  by  Anderson  (26),  Gerster,  et.al.  (27),  Welch  (28),  Iyer  and  Multi  (29). 
This  model  assumes  that  the  rate  of  mixing  or  dispersion  of  a  component  is 
proportional  to  the  concentration  gradient  of  that  component. 

Foss  (30)  has  shown  how  the  residence  time  distribution  function  can  be 
used  to  include  liquid  mixing  effects  in  the  calculation  of  plate  efficiency. 

Mutzenberg  (31)  has  used  the  impulse  test  to  measure  the  residence  time 
distribution  of  a  flow  pattern  on  a  sieve  tray.   He  has  pointed  out  that  the 
total  tray  efficiency  on  a  sieve  plate  based  on  the  vapor  phase  is  a  function 
of  the  local  tray  efficiency,  fluid  velocity,  slope  of  the  equilibrium  curve, 
and  number  of  pools  connected  in  series. 

Leonard  (11)  concluded  that  the  macroscopic  splashing  model  of  Johnson 
and  Maranjosis  (25)  is  unsuitable  for  predicting  dispersion  or  tray  efficiency. 

Olson  and  Heideger  (32)  demonstrated  that  eddy  diffusion  does  not  ade- 
quately describe  the  mixing  on  a  sieve  tray.  He  concluded  that  the  mixing 
is  apparently  of  such  large  scale  as  to  cause  interactions  between  local 
regions  on  the  tray. 

Strand  (33)  measured  the  efficiency  of  desorption  of  acetone,  methyl 
isobutyl  ketone  and  isopropyl  alcohol  from  water  by  air  on  bubble  cap  trays. 
He  derived  an  equation  for  evaluating  the  tray  efficiency.  By  using  an 
axial  dispersed  plug  flow  model  around  which  a  fraction  of  the  entering 
liquid  is  effectively  bypassed,  a  better  agreement  between  experimental  and 
predicted  efficiencies  was  obtained  as  compared  with  the  result  obtained 
from  the  A.  I.  Ch.  E.  method  (34). 

The  above  survey  indicates  that  the  mixing  on  a  distillation  tray  may 
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be  represented  by  either  the  axial  dispersed  model  with  bypassing  or  the 
pool  model  in  series  with  bypassing. 

Cha  (35)  has  discussed  the  ^-distribution  model  at  great  length.   The 
r7 -distribution  model  is  a  mathematical  model  which  is  more  flexible  than 
the  pools  in-series  model.   The  reason  is  discussed  in  the  next  section. 

Johnson  (1)  was  the  first  who  used  the  ^-distribution  model  to  des- 
cribe a  distillation  tray.   Modification  by  addition  of  a  bypassing  stream 
to  the  /'-distribution  was  proved  to  be  satisfactory.   In  this  work,  a 
similar  model,  based  on  previous  work  and  the  physical  reality  of  the  speci- 
fic nonideal  sieve  tray  is  proposed.  An  alternative  model  with  a  different 
path  for  the  bypassing  stream  is  suggested,  which  may  prove  to  be  useful  in 
understanding  the  bypassing  on  the  distillation  tray. 

The  relationship  between  the  /^-distribution  model  and  the  pools  in-series 
model 

The  pools- in-series  model,  or  the  completely  mixed  tanks- in-series 
model  is  an  alternative  approach  to  the  dispersion  model  for  dealing  with 
small  deviations  from  plug  flow. 

Consider  a  vessel  consisting  of  j  pools  in  series.  Then  the  exit  age 
distribution  can  be  obtained  by  a  material  balance  as  (19) 

E(8)  =   J  9   -  e"j9  (6-1) 

(j  -  1)1 

According  to  Cha  (35),  the  mathematical  representation  of  the  /'-distribution 
model  is 

E(t)  =  — i (t  -  D)P"1exp(-  ^A  (6-2a) 
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where 


P(P)  =  /V"1  e-XdX 


it9   p  and  D  are  the  parameters  of  the  model.  D  corresponds  to  the  dead  time 
or  delay  time  of  the  system.  Cha  (35)  also  shown  that 

U   a  V  P  +  D 

and  for  a  closed  system  the  mean  of  E(t)t  denoted  by  u  9   can  be  represented  by 


V 
v 


A    -  t  -  •■ 


Therefore,  one  can  write 


j/i  ^ 


This  means  that  only  two  parameters  among jSf   P,  and  D  are  independent. 
Equation  (6-2)  can  now  be  written  as 

P 

E(t)  =  -= £— (t  -  D)P"1exp[-P(=-^)]  (6-2b) 

(t  -  D)P  f*(p)  L    t  -  D  J 

t  — 

letting  9  =-=r  and  noting  that  E(0)  =  tE(t),  one  obtains 

p  e  -  £• 

E(9)  =  £- (9  .r)P"1exP[-P(r— ^)]  (6-3) 

(1  -Td)   (P)      d  x      rd 


where 


D 
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If  one  considers *r,  as  the  dimensionless  dead  time,  by  letting 'f   =  0 
Equation  (6-3)  reduces  to 

Therefore,  comparing  this  with  Equation  (6-1),  one  can  see  that  the  age  dis- 
tribution function  of  the  tanks- in-series  model  is  identical  to  the  P-dlstri- 
bution  model,  if  the  two  parameters  in  the  P-distribution  model  are  restricted 
by  equating  T*  to  zero,  and  making  p  take  a  positive  integer. 

The  P-distribution  model  with  by-passing  across  the  /^mixing  system  is 
shown  in  Fig.  7.   This  mixed  model  can  be  derived  as  follows: 

Let 

v  V 

T  T 


(6-5) 


and  thus 


v,   v„   Vo 
12    3 

V+V+V    =V  —  +  —  +  ' —  =   1 

Vl  +  V2  +  V3     VT   •       VT   VT   VT 


V2  V3 

v~  s     f2Td   •     T     =  f3rd   •    with  f2  +  f3  =  l    (6"6) 


f~  and  ^3  are  constant. 

Also  let 

V  V  V  V      7(1  -  T  ) 

-             1      1   T  T      K         Cd;  et.   ,. 

t   =j^p  =  —  =  - =  (6-7) 

VI  VT  Vl  VT         3 
and  thus 

ta  -rd) 

V   = —  (6-8) 

0P 


t  is  the  mean  residence  time  based  on  the  entire  system. 
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The  P-mixing  system  without  dead  time  can  be  obtained  by  setting  D  =  0 
in  Equation  (6-2a)  as 

E(t)_,  =  — - tP"1  exp(-  £)  (6-9a) 

Substituting  Equation  (6-8)  into  Equation  (6-9a)  yields 


E(t)   =  — tP"1  exp[_     ]  (6-9b) 


L   BP   J  3P 


Recalling  that  E(0)  ■  tE(t),  one  has 

P  P 

E(9)   =  2-2 QP~l   exp(-  ■  PP   9)  (6-9c) 

n  (l.t//7(P)  l"rd 

Taking  Laplace  transform  of  Equation  (6-9c)  yields 

.      i(s)     =      /""    r L_]p  .  r ££ ]p        (6.9d) 

"    p       (i -rd)pLs+rlp-J       La-rd)s  +  PPJ 

From  the  relation,  E(s)  =  H(s)   ,  one  can  write 

—        C,>  —     —        P  8       P 

E(s)   =  J-  or      C,  =  C,  [t; TTT ]        (6"10) 

-    p      c.  2         l  L  (i  -rd)s  +  p3j 

A  material  balance  over  the  differential  volume  in  the  first  plug  flow 
section  gives 


dC     dC 
"  VTdv     dt 


Rewriting  in  terms  of  reduced  time 
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dC     dC  _t_  -     vx 

Tdv  "  dt  t   *  v^ 

(6-lla) 

v  d£  _  dC 

T  dv  "  d9 


Taking  the  Laplace  transform  and  noting  that  C(0)  =  0  yields 


dC       — 
VT  ~  =  sC  (6-llb) 


Integrating  over  the  boundary  conditions,  C  =  C  at  V  =  0;   C.  =  C.  at  V  =  V,, 
yields 

/*«  =  -f  /2dV 
CQ   C        VT  0 

**T    »  -f  V2  =  -sCf2^d3 
c0        T 

Thus 

Cx  =  C0exp[-  f2srd]  (6-12) 

Material  balance  over  the  second  plug  section,  using  the  same  argument,  gives 
similarly 

c       —  v 

4  dC       s   ,v2 

c3  c     vT  0 

Thus  one  obtains 

C4  =  C3  exp[-  f3srd]  (6-13) 

Material  balance  over  the  second  pot  gives 
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VlC2  +  V2Cl   =  VTC3 


v  v 

1  -  2  - 

C_   =  —  C-  +  —  C. 

T  T 


=  pc2  +  (l  -  3)  ct 


(6-14) 


Substituting  Equations  (6-10),  (6-12),  and  (6-14)  into  Equation  (6-13)  leads 
to  _ 

E(.)  =  Ife)  =  |  =  [pf__£JL_)P  +  (1  .  „,]  exp[-  ,rj     <6-15) 

Taking  the  inverse  transform  yields 

gP+1PP         ro      ,-,'-1        r        PB 


E<9>     =    E— £ [9  -  Td]        exp[-  -i£-  9]  +  (1  -  p)6[9  -  r„]  ; 

>  Td        (6-!6) 


(1  -  TdV     (P)  l  '  Si 


E(0)  =  0  ;     e  <  Td 

Note  that  the  result  is  the  same  as  putting  the  two  equal  sized  plug  flow 
sections  together  as  expected  (1), 

The  /^-distribution  model  with  by-passing  from  inlet  to  outlet  (Fig.  8) 
can  be  derived  in  the  same  way  as  the  first  model.   Let 

v  V 

1  1 

3    =    —      •  1  -Tri    =    — 

H  v  d  V 

T  T 

v,        V0        V„ 
12  3 

12  3  T        »  VT       VT       VT 

V  v 

F    -     f2Td       •       F    =     £3Td  wlth     f2+£3     ■     l 

T  T 
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fj  and  f»  are  constant.  Again,  letting 

\       «"  -  V 

t    =  j/p  =  —  =  

vi  e 

one  has 

t(i  -  r  ) 

t  is  the  mean  residence  time  based  on  the  entire  system. 

The  Laplace  transform  of  the  r-mixing  unit  without  dead  time  is  repre- 
sented by 

b(.)     =     pPg"     [_! ]p  =  [ ES ]p 

_    r      a-rd)pLs+rgpF  La-rd)8  +  ppj 


Since 


one  has 


E(s)   =  J- 
P  C 


C.   =:  C.  T — 1P  (6-17) 


'2 


s  c,  r — 1 

lL(l  -T  )a  +  P3J 


A  material  balance  over  the  first  plug  flow  section  yields,  in  the  Laplace 
transformed  domain, 

f2 
cl  ■  co  exP[-  J  srd  ]  (6"18) 

Similarly  for  the  second  plug  flow  section 

f3 
C3  =  C2  exp[-  jsrj  (6-19) 

Material  balance  over  the  second  pot  gives 
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V1C3  +  V2C0  =  VTC4 

Vl  -    V2  - 
c4  =  -i  c3  +  —  CQ  =  3  c3  +  (I  -  3)  CQ  (6-20) 

T       T 

Substituting  Equations  (6-17),  (6-18),  and  (6-19)  into  Equation  (6-20)  yields 

s  f 

C,   =  Cn  0< LE )P  exp[ 4]   +  (1  -  p)  (6-21) 

°    (I  -Tri)s  +  P3      L   3  J 


Therefore,  the  exit  age  distribution  function  is 

E(s)   =  H(s)   =  ~ 

C0 

sT 

=  3  <t: — z^r )P  exP["  —  ]  +  O  -  P>  C6-22) 

(i  -  r  )s  +  p3         L     3 


Taking  the  inverse  transform 

p+1  p        M 

E(e)   .  _L_E {e.^)P-1«p[-7Jl£_] 

(i  -rd)p/*<p)         p  (1  -Td> 


+  (i  -  3)  6(e)     ?  e   ^  -~ 

3 


E(e)  =  o  ;  e  <  .1 


6  (9)  stands  for  the  delta  function. 


(6-23) 
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CHAPTER  VII 


EXPERIMENTAL 


1.      Equipment 

The  tray  used  in  this  work  is  a  12"  x  48"  rectangular  sieve  tray,  con- 
sisting of  five  zones: 

(1)  Active  or  vapor-dispersion  zone. 

(2)  Distributing  zone. 

(3)  Disengaging  zone. 

(4)  Downcomer  zone. 

(5)  Periphery  waste  zone. 

The  area  occupied  by  the  five  zones  is  shown  in  Fig.  9.   The  spacing  of 
perforations,  which  were  on  an  equilateral-triangular  basis,  range  from  2.5  to 
5  diameters.   With  spacing  closer  than  2.5  diameters,  excessive  coalescence  of 
bubbles  occurred.  At  spacings  greater  than  five  diameters,  tray  capacity  was 
limited  and  inactive  liquid  zones  developed  between  the  perforations,  which 
probably  reduced  transfer  efficiency.   The  percentage  of  free  area,  which  is 
defined  as  the  ratio  of  the  perforation  area  to  the  active  area,  was  designed 
to  be  far  less  than  the  commercial  range  of  5  to  15%. 

The  gauges  for  measuring  the  liquid  depth  and  the  static  pressure  between 
the  plates  was  specially  designed  as  shown  in  Fig.  10. 

The  accessories  and  the  overall  set  up 

The  overall  experimental  set-up  is  shown  in  Fig.  11.   The  constant  head 
tank  shown  in  the  upper  left  hand  corner  was  a  150  gallon,  stainless  steel, 
cone-bottom  vessel.   The  tank  was  supplied  with  tap  water  from  the  building 
main  at  a  rate  of  30  gallons  per  minute  and  provided  a  constant  head  of  20  feet. 
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Two  Brooks  Fullview  rotameters  were  used  to  meter  the  water  to  the  test 
trays.   These  rotameters  had  ranges  of  0.5  to  5.4  g.p.m.  and  2.0  to  23  g.p.m. 
The  guaranteed  accuracy  of  these  meters  was  ±  2%  of  maximum  scale  reading  from 
10%  to  100%  of  flow,  but  the  check  calibrations  indicated  that  flow  rates 
could  be  read  accurately  to  within  1.5%  of  the  scale  reading.   The  calibra- 
tion curves  are  given  in  the  Appendix  la  and  lb. 

One  Brooks  Fullview  rotameter  with  a  range  of  5  -  50  s.c.f.m.  was  used 
to  meter  the  air  to  the  test  tray.   The  calibration  of  the  air  rotameter  was 
checked  with  a  wet  test  meter.   The  temperature  and  pressure  calibration 
correction  curves  are  given  in  Appendix  2a,  2b. 

In  order  to  get  a  uniform  inlet  concentration  of  the  tracer  through  the 
cross-section,  a  distributor  was  provided  so  that  the  water  was  divided  into 
five  streams  which  were  evenly  distributed.   The  five  streams  flowed  into  the 
downcomer  which  was  packed  with  spherical  glass  particles.   The  glass  particles 
were  used  to  obtain  uniform  flow  and  to  reduce  bubbling  of  the  entering  fluid 
(see  Fig.  12). 

High  pressure  air  was  fed  into  a  30  gallon  surge  tank  where  it  was  bubbled 
through  water  in  order  to  saturate  it.   The  combination  of  the  pressure  regu- 
lator and  the  surge  tank  served  to  damp  out  noticeable  fluctuations  in  air 
pressure. 

2.   Instrument  of  measurement 

The  electrical  conductivity  method  was  employed  in  this  work.   Continuous 
recording  of  conductivity  has  been  widely  used  by  many  previous  researchers 
(39,  40,  41,  42).   Prausnitz  (43)  used  a  conductivity  probe  to  measure  the 
turbulent  mixing  properties  and  Manning  (44)  improved  the  conductivity  probe 
to  measure  the  concentration  fluctuations  in  a  baffled  stirred  vessel. 
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Mutzenberg  (31)  used  the  conductivity  instrument  to  measure  the  residence 
time  distribution  on  a  sieve  tray. 

The  conductivity  measurement 

Electrolytic  conductivity  is  a  measure  of  the  ability  of  a  solution  to 
carry  an  electric  current.   The  specific  conductance  is  defined  as  the  reci- 
procal of  the  resistance  in  ohms  of  a  1  cm  cube  of  the  liquid  at  a  specified 
temperature.   The  commonly  used  units  of  specific  conductance  are  the  reci- 
procal ohm-cm  (or  mho/cm),  and  one  millionth  of  this,  micro-mho/cm.   High 
quality  condensed  steam  and  distilled  or  deminerallized  water  have  specific 
conductances  at  room  temperature  as  low  as,  or  lower  than,  1  micro-mho/cm. 
This  would  be  caused  by  about  one-half  ppm  of  a  dissolved  salt.   (sodium 
chloride  solution  was  used  as  a  tracer  in  this  work).  Absolutely  pure  water 
has  a  conductivity  of  0.055  micromhos/cm  at  25°C. 

The  mechanism  of  electrolytic  conduction  differs  from  that  of  metals. 
In  the  latter,  current  is  carried  by  free  electrons;  in  the  former,  current 
is  carried  by  ions  only.   Electrolytic  conductivity  is  non-specific  and  all 
ions  present  in  a  solution  contribute.   This  is  in  marked  contrast  to  the 
general  lack  of  influence  of  other  ions  on  pH  measurement  which,  over  most 
of  its  useful  range,  is  responsive  to  the  H-ion  concentration  alone. 

The  conductivity  probe 

In  this  work  small  conductivity  probes  were  desirable  in  order  to  mini- 
mize disturbance  of  the  flow  pattern  by  the  probes.   Two  probes  designed  by 
the  Industrial  Instrument  Inc.  were  used  (Appendix  3).   These  two  probes 
consisted  of  two  small  metal  plates  situated  at  the  end  of  a  5  ro.m.  outside 
diameter,  12"  long  epoxy  tube.  The  exposed  ends  of  these  two  plates  of 
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metal  were  platinized  (Fig.  13). 

The  conductivity  meter 

The  model  RA4  Solu  Meter  used  is  a  direct  reading  instrument  providing 
a  continuous  indication  of  solution  conductivity.  A  0  to  10  DC  signal  pro- 
portional to  the  conductivity  is  provided  to  operate  a  remote  potential 
recorder.  A  non-bridge  type  of  measuring  circuit  is  utilized.   The  conduct!" 
vity  ranges  from  0  to  2500  micromhos/cm.   The  specifications  of  this  conduc- 
tivity meter  are  shown  in  Table  1.   Fig.  14a  shows  the  block  diagram  in  the 
instrumental  set  up.   Fig.  14b  is  a  photograph  of  the  instrumental  set  up. 
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To    solu   meter 


5  m.  m.dtGinofar 

i  2"  long    opoay   tuba 


Tw/o    plotiniiod      plofos 


Fig.  13.  The     conductivity     probe. 
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TABLE   1.      SPECIFICATIONS  OF  TYPE   RA4  SOLU  METER 


Ranges 

Temperature  compensation 
Cell  voltage 
Measurement 

Output 

Accuracy  of  output 

Power  requirements 

Mount  ing 

Case 

Dimensions 

Weight 


Three  standards  2000,  200  and  100  ohms  measured 

resistance  for  full  scale  deflection* 

Mannual  or  automatic 

3  volts  maximum 

In  undergrounded  or  grounded  electrolytic 

solutions. 

0-10  mv  DC  across  100  ohms 

±1%  of  full  scale,  above  1/10  full  scale 

18  watts  at  120  volts 

Panel  or  portable 

Aluminum  gray  finish 

9M  high  (over  handle)  x  6-1/4"  wide  x  8"  deep 

8-1/2  pounds 


49 


BecKman    offner 

type  RC  dyno- 
graph  recorder 

Solu    meter 

Conduc- 
tivity 
probe 

Type  RA4-WB-SI7-TI 

Ser.  No.  74  4 

Cell     constant     0.500 /cm. 

Fig.  14.  Block     diagram     of    the   instrumental 


set    up. 
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CHAPTER  VIII 

PHOTOGRAPHIC  STUDY  AND  VISUAL  OBSERVATION 

The  purposes  of  the  photographic  study  and  visual  observation  were  as 
follows: 

(1)  To  provide  a  physical  picture  of  the  flow  pattern  occuring  on  the 
sieve  tray. 

(2)  To  observe  the  effect  of  flow  patterns  (oscillating,  normal  flow 
patterns,  etc.)  on  the  degree  of  liquid  mixing.   The  oscillating  flow  pattern 
of  the  first  type  may  be  especially  important  in  this  respect  since  it  may 
decrease  the  backmixing  and  increase  the  Murphree  vapor  efficiency  (51). 

(3)  To  establish  a  plausible  mechanism  for  a  bypassing  mechanism  on 
the  distillation  tray. 

Formation  of  the  stratified  elements 

When  the  air  flowed  at  a  very  low  rate  (say  2  g.p.m.),  the  frontal  edge 
of  the  potassium  permanganate  solution  injected  as  a  tracer  proceeded  very 
uniformly  near  the  inlet.  (It  is  convenient  to  use  KMnO,  solution  as  a  tracer 
for  photographic  study,  since  it  has  color).  At  a  position  about  half  way 
between  the  inlet  and  the  bubble  or  active  zone,  finger like  projections  of 
the  dye  front  occurred  predominantly  at  or  near  the  center.   This  fingerlike 
projections  (probably  stratified  elements)  appeared  to  be  drawn  forward 
through  the  action  of  the  bubble  fountain,  producing  a  non-uniform  dye  con- 
centration throughout  the  cross-section  (Fig.  15).  At  a  position  half  way 
between  the  inlet  and  the  bubble  zone,  the  concentration  of  the  dye  tracer 
appeared  to  be  smaller  in  the  vicinity  of  the  walls,  where  the  velocity  of 
the  stream  flow  was  nearly  equal  to  zero.   The  formation  of  the  stratified 
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Stage  one 


Stage  two 


1 


Fig.  15  a, 


Bubble  zone  Tracer  flows  in 

or 

Fountain  zone  uniformly 


Direction  of 
flow        ^_ 


Fig.  15  b. 


Occurrence    of  non-uniform 
distribution    of  tracer. 


Stage  three 


Fig.  15  a 


Tracer  flows    out   nearly  uniformly 
Fig.  15.  Stages    in  the  formation   of    stratified   elements 
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elements  might  be  explained  by  on  increase  in  the  non-uniform  velocity 
profile  which  in  turn  enhanced  by  the  non-uniformity  in  tracer  distribution 
which  already  existed  at  the  inlet. 

The  frontal  edge  of  the  dye  tracer  at  the  bubble  zone  and  at  the  outlet 
again  appeared  to  be  uniform,  apparently  due  to  the  stirring  and  agitation 
of  the  air  bubbles. 

Figs.  16a,  16b,  16c,  and  16d  show  the  formation  of  stratified  elements 
under  various  conditions. 

Formation  of  cellular  foam  and  bubble  bursting 

Here  observations  were  confined  to  operations  with  a  multiply  perforated 
tray.  The  hole  diameter  was  3/16"  and  the  spacing  was  1-1/2".  Fig.  17.  is 
a  photograph  taken  under  the  conditions  with  the  air  flow  rate  of  7.5  c.f.m., 
water  flow  rate  of  7.5  g.p.m. ,  weir  height  of  3  in.  and  7  rows  of  holes.  One 
can  see  the  formation  of  bubbles,  the  growth  of  bubbles,  and  the  bursting 
of  bubbles. 

When  the  air  velocity  was  less  than  approximately  10  c.f.m.,  the  fre- 
quency of  the  bubbling  was  low  and  hence  sufficient  time  became  available 
for  the  growth  of  the  bubbles.   When  a  bubble  collided  with  others,  two 
events  might  happen:  one  is  the  formation  of  a  larger  bubble  and  the  other 
a  breakdown  into  smaller  bubbles.  At  the  low  air  velocity,  the  formation  of 
the  larger  bubble  would  be  more  probable,  especially  near  the  surface.   Hence 
a  collection  of  the  larger  bubbles  flows  upwards.  When  it  reaches  the  sur- 
face, it  becomes  cellular  foam  (Fig.  18).  Of  course,  the  formation  of  the 
cellular  foam  not  only  depended  on  the  air  velocity,  but  also  on  the  liquid 
properties,  such  as  surface  tension.  At  a  very  high  air  velocity (approxi- 
mately over  20  c.f.m.  in  this  system),  the  frequency  of  the  bubbling  was 
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quite  high.  Sufficient  time  was  probably  not  available  for  the  growth  of 
the  bubbles  at  high  air  flow  rates.  The  bubbles  at  a  high  air  flow  rates 
were  easily  broken,  and  a  collection  of  smaller  bubbles  was  formed  at  the 
top  of  the  bubble  zone,  which  appeared  as  froth. 

Studies  of  the  relations  between  the  bubble  formation  and  mass  transfer 
can  be  found  in  the  published  literature  (46,47,48), 

Formation  of  the  fountain  and  bypassing  effect 

The  behavior  of  the  fountain  or  bubble  zone  may  be  qualitatively  ex- 
plained as  follows.   In  a  circular  jet,  the  volume  of  flow  at  a  given  dis- 
tance from  the  orifice  is  independent  of  the  momentum  or  the  excess  of  pressure 
under  which  the  jet  leaves  the  orifice.  A  jet  which  leaves  under  a  large 
pressure  difference  remains  narrower  than  one  leaving  with  a  small  pressure 
difference.  The  latter  carries  with  it  comparatively  more  stationary  fluid 
than  the  former. 

The  stationary  fluid,  after  the  injection  of  air,  no  longer  remains 
stationary.  The  horizontal  velocity  (or  velocity  in  the  radial  direction 
perpendicular  to  the  flow  through  the  jet)  of  water  is  increased  at  a  dis- 
tance near  the  bottom  of  the  wall;  while  a  part  of  the  stationary  fluid,  far 
from  the  wall  acquires  both  horizontal  and  vertical  (or  axial  direction,  the 
same  direction  as  flow  through  a  jet)  velocities.   If  the  stationary  fluid 
(water)  has  its  boundary  exposed  to  the  air,  part  of  the  fluid  (water)  will 
possibly  be  carried  to  the  top  of  the  jet*   Figure  19a  shows  the  change  in 
velocity  profile  with  respect  to  distance  and  the  stream  line  pattern  of  a 
circular  jet. 

The  behavior  of  the  flow  pattern  on  the  sieve  tray  may  be  explained  by 
means  of  an  analogy  to  a  circular  jet.  The  injected  fluid  was  air  and  the 
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Fig.   19a    Stream    line    pattern    for    a    circular  lamilar  jet(45). 
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interface 
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Fig.  19b.    Hypothetical   flow    pattern    for  air  passing  through 
a   distillation   tray    with  one- row  of    holes. 
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liquid  on  the  tray  was  actually  not  stationary.   However,  the  holes  of  the 
tray  may  be  regarded  as  a  group  of  jets.   Figure  19b  shows  the  hypothetical 
flow  pattern  for  air  passing  through  a  distillation  tray  with  one  row  of 
holes  perpendicular  to  the  direction  of  flow.   When  the  upwards  liquid  flow 
reaches  the  interface  (interface  between  the  vapor  space  and  liquid),  it  is 
then  divided  into  two  streams,  one  stream  flowing  to  the  direction  of  the 
original  liquid  flow,  and  the  other  flowing  against  the  direction  of  the 
liquid  flow  as  shown  in  Fig.  19b.   This  gives  rise  to  the  formation  of 
fountain.   The  phenomenon  of  the  fountains  will  be  discussed  later  in  con- 
nection with  the  swinging  flow  pattern  or  oscillating  flow  pattern  of  the 
second  type. 

The  mechanism  for  the  occurrence  of  bypassing  streams  on  a  distillation 
tray  still  remains  unknown.  Mutzenberg  (31)  carried  out  an  interesting 
investigation  in  which  he  compared  the  tray  efficiencies  evaluated  under 
different  conditions.   He  found  that  the  tray  efficiency  was  increased  by 
20%  when  baffle  plates  were  inserted  into  the  bubbling  zone  at  different 
depths.   This  may  enable  one  to  suspect  that  bypassing  may  occur  through  the 
froth  and  can  be  reduced  by  inserting  the  plates.  Another  explanation  is 
that  the  residence  time  of  the  fluid  element  is  increased  when  several  plates 
are  inserted  into  the  bubble  zone.   If  one  fit  the  data  both  with  and  without 
the  plates  to  the  F-distribution  model  with  bypassing  across  the  f-mixing 
unit,  the  result  might  show  a  change  in  the  bypassing  parameter,  and  hence 
indicate  which  is  the  more  plausible  mechanism. 

Occurrence  of  oscillating  flow  pattern 

The  four  variables  investigated  in  this  work  were  weir  height,  liquid 
flow  rate,  gas  flow  rate,  and  the  percentage  free  area  for  gas  flow. 
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The  oscillating  flow  patterns  observed  could  be  classified  into  two 
types,  one  that  oscillates  perpendicular  to  the  direction  of  liquid  flow 
(the  first  type)  and  the  other  parallel  to  the  direction  of  liquid  flow 
(the  second  type). 

In  runs  where  the  weir  height  was  1",  normal  flow  patterns  were  obtained 
regardless  of  the  changes  made  in  the  other  variables  within  the  range  of 
interest  (gas  flow  rate  of  approximately  5  to  20  c.f.m.,  liquid  rate  of 
approximately  5  to  15  g.p.m. ).   The  term,  normal  flow  pattern,  is  used  here 
to  distinguish  it  from  an  oscillating  flow  patterns. 

However,  at  a  weir  height  of  1.5  in.  and  a  liquid  rate  of  10  g.p.m.  with 
a  gas  rate  of  10  c.f.m.  and  3  rows  of  holes,  an  oscillating  phenomenon  was 
observed.   Careful  observation  showed  that  the  flow  pattern  fluctuated  regu- 
larly.  The  peaks  rose  up  the  side  walls  leaving  a  central  valley  and  then 
damped  down.  As  this  happened  a  peak  was  formed  in  the  approximately  central 
region.   The  process  was  repeated  regularly.   When  the  peaks  on  both  sides 
were  decreasing  in  height,  water  forced  down  through  the  holes  near  the  sides: 
while  a  peak  in  the  center  was  diminishing,  water  came  through  holes  near  the 
sides.   The  amplitude  of  the  oscillation  was  larger  at  higher  air  flow  rates 
and  liquid  flow  rates.   Figure  20  illustrates  the  stages  of  the  first  type 
of  oscillating  flow  pattern.   While  the  oscillating  flow  pattern  shown  in 
Fig.  21a, b  was  slightly  different  from  that  depicted  in  Fig.  20,  they  essen- 
tially show  the  oscillating  flow  pattern  of  the  first  type.   The  front  view 
of  the  normal  flow  pattern  is  given  in  Figure  20c  for  comparison.   Figures  21a 
and  21b  were  taken  under  the  same  conditions.   Figures  22a, b,c,  and  d  show 
the  side  view  of  the  oscillating  flow  pattern  of  the  first  type. 

Figures  23a  and  b  show  the  side  view  of  the  swinging  flow  pattern,  that 
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Fig.  20.  Stages    in   the  formation   of 
ocillating    flow    pattern, 
(front  view,  facing  in  tire  direction 
from   downstream   to  upstream). 
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is,  the  second  type  of  oscillating  flow  patterns.   The  swinging  flow  pattern 
occurred  at  very  low  liquid  flow  rate  (2  g.p.m.  or  less).   The  direction  of 
the  oscillation  was  different  from  the  first  type  oscillating  flow  pattern. 
It  oscillated  parallel  to  the  direction  of  the  liquid  flow.   The  swinging 
type  oscillation  phenomenon  could  be  eliminated  by  simply  increasing  the 
liquid  flow  rate.   The  side  view  of  the  normal  flow  pattern  is  given  in  Fig. 
23c  for  comparison. 
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Explanation  for  the  occurrence  of  the  oscillating  flow  patterns 

A  possible  explanation  for  the  occurrence  of  the  oscillating  flow  pat- 
tern of  the  first  type  is  given  below: 

Due  to  the  non-uniform  velocity  profile  (the  water  flowing  near  the 
center  region  of  a  tray  being  faster  than  that  near  the  wall),  a  static  pres- 
sure difference  may  exist  between  the  central  region  and  the  region  near  the 
wall.  When  the  upwards  flow  of  air  is  just  at  its  critical  pressure  to 
overcome  the  overall  liquid  head,  it  tends  to  pass  through  the  holes  where 
the  pressure  head  is  low.   Supposing  that  the  low  pressure  head  first  occurs 
in  the  central  region  where  the  water  flow  velocity  is  higher  than  that  at 
both  sides  as  mentioned  above,  the  air  stream  will  cause  the  liquid  level 
to  peak  in  the  center  region.  As  the  peak  reaches  its  maximum  height,  the 
region  of  the  low  pressure  head  moves  toward  one  or  both  sides  near  the  wall. 
Now  the  air  tends  to  flow  through  the  region  near  the  wall  and  this  in  turn 
causes  a  peak  or  peaks  to  be  formed  in  the  corresponding  region  as  illustrated 
in  Fig.  20.   The  repeated  processes  cause  the  oscillation  phenomenon. 

The  second  type  of  oscillating  flow  pattern  might  be  explained  as  below: 
High  velocity  air  (air  flow  rate  of  approximately  ten  c.f.m.  for  from 
one  to  three  rows  of  holes  in  the  sieve  tray  under  consideration)  can  readily 
overcome  the  liquid  depth.   The  flow  of  liquid  under  such  a  high  air  flow 
rate  tends  to  be  blocked  by  the  bubble  zone.   This  in  turn  tends  to  cause 
the  liquid  to  accumulate  temporarily  behind  the  bubble  zone  until  its  level 
reaches  a  certain  depth  to  overcome  the  "blocking"  action  by  the  bubble  zone. 
The  repetition  of  this  process  gives  rise  to  the  oscillation  phenomenon  of 
the  second  type. 
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CHAPTER  IX 

QUANTITATIVE  ANALYSIS 

The  response  of  the  measuring  instrument  must  be  understood  before  the 
details  of  the  quantitative  analysis  of  the  experimental  data  can  be  carried 
out. 

Range  of  the  conductivity 

The  range  of  the  conductivity  of  common  solutions  is  shown  in  Fig.  24. 
The  minimum  concentration  detected  in  this  work  is  the  concentration  of  tap 
water.   The  salt  solution  being  measured  here  was  a  function  of  the  water 
flow  rate,  and  the  concentration  and  amount  of  the  salt  solution  which  was 
injected.   If  the  range  were  so  chosen  that  it  was  much  over  the  required 
range,  the  output  signal  of  the  Solu  Meter  would  be  very  weak  so  that  the 
recorder  would  fail  to  strengthen  it.   On  the  other  hand,  if  the  range  was 
too  small  the  variation  of  the  tap  water  salt  concentration  would  have  intro- 
duced appreciable  error  into  the  recorded  curves.   For  this  reason,  the  daily 
variation  of  the  tap  water  salt  conductivity  was  measured  continuously  for 
one  month.   There  was  no  significant  change. 

Problems  of  the  linearity  of  the  relation  between  concentration  and 
conductivity  and  temperature  effect  on  the  conductivity  change  have  been 
partially  discussed  in  Chapter  7.   These  two  problems  are  discussed  further 
below. 

Relationship  between  conductivity  and  concentration 

The  relationship  between  conductivity  and  concentration  is,  in  general, 
non-linear.   However,  at  very  low  concentration  it  is  approximately  linear. 
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Prausnitz  (43)  has  presented  the  following  equation  for  his  specific  instrument* 

8.73  -  0.136c  (9.u 

R0.93 

where 

N  =  normality  of  the  hydrochloric  acid, 

R  =  resistance,  ohm, 

t  =  temperature,  °C. 

In  this  experiment  the  absolute  concentration  was  not  required.   Two 
concentration-time  responses,  one  at  the  inlet  and  the  other  at  the  outlet, 
were  needed.   The  important  requirements  were  that 

(1)  The  recorded  millivolts  be  proportional  to  the  conductivity.   It 
was  necessary  that  the  conductivity  be  proportional  to  the  concentration 
within  the  range  of  interest. 

(2)  The  two  probes  have  the  same  efficiencies. 

The  first  requirement  was  of  the  utmost  importance,  while  the  second 
was  not  absolutely  necessary  since  a  correcting  program  could  adjust  the  data. 
Therefore  it  was  necessary  to  make  two  equipment  tests,  one  before  and  the 
other  after  each  experimental  run.   This  primary  test  was  carried  out  as 
follows: 

(1)  Prepare  1000  c.c.  of  the  tap  water  in  a  container  and  immerse  the 
probes  into  the  water.  Adjust  the  recorder  pens  to  get  the  base  zero  line 
(The  tap  water  is  collected  directly  from  the  system). 

(2)  Drop  one  c.c.  of  the  salt  solution  into  the  container  and  stir  it 
thoroughly. 

(3)  Observe  the  recorded  curve  until  it  reaches  the  plateau.   Then 
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record  the  change  in  the  millivolt  reading. 

(4)  Drop  another  c.c.  of  the  same  salt  solution  and  record  the  change 
in  the  millivolt  reading. 

(5)  Repeat  the  procedure  until  the  range  of  interest  has  been  exceeded. 
The  results  first  obtained  for  one  of  the  probes  were  not  satisfactory, 

but  after  careful  adjustment  of  the  manual  temperature  compensation,  two 
linear  curves  were  obtained.  The  test  was  repeated  after  each  run  to  make 
sure  that  the  temperature  change  during  the  period  of  operation  did  not  change 
the  recording.  Two  typical  curves  are  given  in  Figs.  25  and  26. 

Since  the  amount  of  salt  solution  introduced  was  not  over  10  c.c.,  the 
concentration  increase  could  be  treated  as  an  arithmetical  progression. 

Experimental  data 

Table  2  shows  a  sample  of  the  data  sheets  used  in  recording  the  experi- 
mental data.  Two  sets  of  data  were  taken  in  order  to  analyze  the  effect  of 
the  stratified  elements  on  the  degree  of  mixing.  One  set  of  data  was  obtained 
using  one  row  of  holes  (15  holes)  in  the  bubble  zone  while  the  other  set  was 
collected  using  three  rows  of  holes  (43  holes).   Figs.  27a  and  27b  show 
records  of  the  input  and  output  pulses  obtained  in  two  duplicates  of  run  11. 
Both  the  input  and  the  output  concentration- time  curves  appear  to  be  very 
similar.   This  gives  some  assurance  that  the  recording  devices  were  function- 
ing normally. 

Figures  28a,  28b,  and  28c  show  records  of  the  input  and  output  pulses 
for  three  runs  (runs  12,  19,  and  20)  at  different  water  flow  rates  but  at  a 
constant  air  flow  rate  of  8  c.f.m.  Appearance  of  peaks  on  the  response  curves 
for  these  runs  are  appreciably  different.   Run  12  was  carried  out  at  a  water 
rate  of  10  g.p.m.  and  an  air  rate  of  8  c.f.m.   Run  19  was  carried  out  at  a 
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Fig.  25 
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Primary    test   for   the    linearity 
of  the    infet    probe   before 
experiment  runs. 
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Fig.  25b.  Primary  test  for  the  linearity 
of  the  outlet  probe  before 
experimental   runs. 
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Fig. 26a.  Primary  test  for  the  linearity  of  the 
inlet  probe  after  experimental   runs. 
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Primary  test  for  the  linearity  of  the 
outlet  probe  after  experimental    runs. 


TABLE   2 


Data  Sheet 
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Run  Number 


11 


Tray  Variables 

Number  of   Holes  15    (1   row) 


Weir  Height 


1.5 


in 


Liquid  Depth    (Average)     2.1292    in 
Pressure  Drop       9.1  cm  H2° 


Date        May   19,    1966 
Water  Flow  Rate        10  g.p.m. 
Air  Flow  Rate  8  c.f.m. 


Water  Temperature  15.5°C 
Room  Temperature  2 5 . 5  °C 
Total  Pressure  Drop  9.1  cm  H2° 


Position  of  the  Probes  Located 


Inlet 


Center 


Outlet   2.4 


in 


from  left 


Local 

Liquid  Depth 

1. 

5.25  cm 

2. 

5.31 

3. 

5.45 

4. 

5.37 

5. 

5.50 

6. 

5.57 

Instrument  Settings 
Scale   • Ix  to  • lx 


Time  Constant    0.5  sec 


Preamplifier     .5/. 5 
Chart  Speed    5  ■*/■•; 


Preamplifier  Gain 


Comments 


Tracer  Volume   10  c.c.  saturated  salt  solution 
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Fig.   27a.      The   input  and  output  curves  of  Run  11. 

27b.     The   input  and  output  curves  of  a  duplicate  run  for  Run  11. 
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Fig.  28a.  The  input  and  output  curves  of  Run  12. 
28b,  The  input  and  output  curves  of  Run  19. 
28c.   The  input  and  output  curves  of  Run  20. 
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very  high  flow  rate  of  water  (16  g.p.m.),  therefore,  the  broadening  of  the 
response  curve  as  compared  with  the  inlet  signal  was  almost  indistinguishable. 
Run  20,  carried  out  at  a  very  low  water  flow  rate  (4  g.p.m.),  gave  a  broader 
peak  with  a  very  long  tail.  The  fluctuations  in  the  curves  in  Run  20  were 
probably  due  to  the  oscillating  flow  patterns  of  the  second  type  (swinging 
back  and  forth  in  the  direction  of  flow)  described  in  the  preceding  chapter. 

As  previously  mentioned,  the  small  uncertainty  as  to  the  end  of  the 
response  curve  causes  a  significant  error  in  the  calculation  of  the  variance 
of  the  response  curve.   This  was  tested  by  sightly  changing  the  cut-off  point 
of  the  tail,  which  caused  a  change  in  the  total  area  of  less  than  1%.   How- 
ever, it  changed  the  variance  about  10%. 

In  addition,  appreciable  differences  between  the  total  areas  of  the 
input  signal  and  output  signal  were  found  in  some  of  the  experimental  runs. 
This  was  probably  partially  due  to  uneven  distribution  of  the  inlet  tracer 
and  also  to  the  stratification  effect.  Therefore*,  a  systematic  adjustment 
of  the  data  was  required. 

Roose  (49)  proposed  the  method  which  changes  the  ordinates  of  Y(9)  with 
the  empirical  correction  formula, 

Y(9)c  =  PY(9)  [1  -  Q9]  ,  (9-2) 

where  Y(9)  represents  the  adjusted  response  data,  Y(9)  the  raw  normalized 
response  data,  P  and  Q  are  the  correction  factors.   If  the  data  are  good,  the 
values  of  P  and  Q  should  be  close  to  one  and  zero  respectively.   For  the 
pulse  testing,  it  can  be  assumed  that  the  input,  X(9),  is  fairly  correct  and 
that  it  is  only  necessary  to  correct  the  output  data.  P  and  Q  can  be  com- 
puted by  noting  two  properties  of  the  pulse  testing  data. 
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The   first   property   is 

CO  oo 

/  x(e)de  =  /  Y(e)de  =  1  (9-3) 

0         0 

The  second  property,  which  has  been  stated,  is  that  the  mean,yWA  ,  of  Y(8) 
minus  the  mean,X/   ,  of  X(9)  is  equal  to  the  mean.//  ,  of  E(9),  that  is, 

S&Y  "/ex  =  V0YX  =//e  (9-4) 

Zwietering  (11)  has  shown  that  for  a  closed  system  the  mean  of  E(9)  is  equal 
to  unity.   Therefore 

/"evVex   =Ve«   "    '  (9"5> 

From  equations  (9-2)  and  (9-3)  we  have 


/  Y(9)  d9  =P/  Y(9)d9  -  PQ  /  8Y(9)d9  =  1  (9-6) 

0      c       0  0 

From  equations  (9-2),  (9-5)  and  (9-6) 


/  9Y(9)  d9  -  /  9X(9)d9 
0       c     0 

OO  OO  f.  oo 

=  P/  9Y(9)d9  -  PQ  /  9  Y(9)  d9  -  /  9X(9)d9  =  1 


By  letting 


oo 

A  =  /  9Y(9)d9  , 
0 


B  =  /  9X(9)d9  ,  (9-8) 

0 


and 
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°°  2 
M  =   /  0  Y(6)de   , 

0 
the  solution  of  equations  (9-6)  and  (9-7)  for  P  and  Q  can  be  obtained  as 


P  -  A(l  +  B)  -  M 

2  * 

A   -  M 

(9-9) 

1   +  B  -  A 


Q  = 


A(l  +  B)  -  M 


Tables  3  and  4  contain  the  conditions  of  the  experimental  runs,  Table  5  and  6 
list  values  of  the  correction  factors  P  and  Q. 

Analysis  of  the  data  and  results 

The  following  procedures  were  used  to  obtain  the  mean  and  variance  from 
the  input  pulse  of  X(t),  and  output  pulse  of  Y(t). 

(1)  Read  the  input  and  output  pulses* 

Since  the  recorded  pulses  were  sufficiently  smooth,  the  time  increment  for 
X(t)  of  0.5  sec.  and  that  for  Y(t)  of  1*0  sec.  were  used.   The  average  points 
counted  for  X(t)  and  Y(t)  were  nearly  the  same.   It  was  about  80  points  res- 
pectively.  Much  smaller  increments  were  necessary  if  the  pulse  was  not 
sufficiently  smooth. 

(2)  Punch  the  data  read  from  the  X(t)  and  Y(t)  curves  on  the  IBM  cards. 

(3)  Calculate  P  and  Q  and  obtain  a  set  of  corrected  data  from  the  data 

correction  program. 

2 

(4)  Use  the  corrected  data  to  calculate  the  ^/^Ayvt  an<*Aa"*AYX  accord" 

ing  to  the  following  equations: 

The  first  moment,  or  the  mean,  of  a  distribution,  f (t),  is  defined  as  (8) 
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TABLE    5.      COK  :    CT.ION   FACTORS    P  ANu  Q    FROM   RUN    1    TO    RUN   8. 


Run  Number 

P 

Q 

1 

2.03129 

0.320165 

2 

2.45992 

0.337048 

3 

1.74853 

0.257272 

4 

1.93885 

0.287961 

5 

2.07940 

0.296839 

6 

2.06114 

0.284637 

7 

1.97809 

0.284637 

8 

2.09876 

0.282301 

average 

2.056 

0.294 

TABLE   6.      CORRECTION  FACTORS   P  AND  Q   FROM  RUN   11   TO   RUN   1'8. 


Run  Number 

P 

Q 

11 

3.53030 

0.382755 

12 

3.19750 

0.380024 

13 

2.85484 

0.377417 

14 

2.34872 

0.358905 

15 

2.21319 

0.326604 

16 

2.93831 

0.363699 

17 

2.91013 

0.369718 

18 

2.78836 

0.368090 

average 

2.847 

0.366 
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/   tf(t)dt 
u     =  2 (9.10) 

/   f(t)dt 
0 

The  seeond  moment  about  the  mean,  or  the  variance,  is  defined  as  (8) 


/  (t  -ii  )  f(t)dt 

2     0      '  c 

a       =  (9-li) 

t        °° 

/   f(t)dt 

0 


If  f (t)  is  read  at  N  equidistant  points,  equations  (9-10)  and  (9-11)  may  be 
expressed  for  the  purposes  of  numerical  computation  as 

N 

Mt  m  ~H (9-12) 

ifif<ti) 


S  (t.  -u)   f(t) 
2     i=l    i    t 

a     = (9-13) 

t         N 

S  f(t.) 
i=l    J 


A  simple  computer  program  was  developed  by  Johnson  (1)  to  compute  the  dimen- 

2 
sionless  quantities,  ^/^qyv  and  *  afiYX#  Tnis  was  done  by  first  computing^/ 

2 
and  a  for  both  X(t)  and  Y(t)  according  to  equations  (8-12)  and  (8-13).  The 

difference  in  the  means  was  then  calculated  and  made  dimensionless  by  dividing 

it  by  t,  the  mean  residence  time.  The  difference  in  the  variances  was  cal- 

-2 
culated  and  then  normalized  by  dividing  it  by  t  ,  The  results  can  be  found 

in  Tables  7,  and  8. 


93 


TABLE   7.      EXPERIMENTAL   VALUES   OFyeyx  AND  a0"         FROM  ADJUSTED 
DATA    FOR   RUN   1    TO    RUN   8. 


Run  Number 

Vgyx  Ve 

2              2 
Aa        =  a 

9YX         e 

1 

1.0009 

.0792 

2 

1.0009 

.0669 

3 

1.0007 

.2052 

4 

1.0007 

.2075 

5 

1.0008 

.1418 

6 

1.0008 

.1608 

7 

1.0007 

.1587 

8 

1.0008 

.1233 

TABLE   8.      EXPERIMENTAL  VALUES  OF^  AND*<7QYX   FR0M  ADJUSTED 

DATA   FOR  RUN   11   TO   RUN   18. 


Run  Number 

Vgyx  Ve 

2              2 
aq        =  a 

0yx       e 

11 

1.0278 

-.1843 

12 

1.0277 

-.0600 

13 

1.0276 

-.0348 

14 

1.0275 

-.0646 

15 

1.0274 

.0906 

16 

1.0276 

-.0043 

17 

1.0276 

-.528 

18 

1.0276 

.0161 
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The  third  moment  about  the  mean,  or  the  skewness,  if  required,  can  be 
calculated  from  its  definition, 


3 

/  (t  -  u  )    f(t)dt 

3     0       c 

t;  = (9-i4) 

t  oo 


/  f(t)dt 
0 


If  f (t)  is  read  at  N  equidistance  points,  equation  (8-14)  may  be  expressed 
for  the  purposes  of  numerical  computation  as 


X  (ti-//t)3f(ti> 


«-  3  =  — ■  (9-15) 

*-  t         N 

£  f(t.) 
i=l    * 

3 
This  can  be  done  by  first  computing  f   for  both  X(t)  and  Y(t)  according  to 

equation  (9-15),  The  difference  in  the  skewness  Can  then  be  calculated  and 

_3 
made  dimensionless  by  dividing  it  by  t  . 


Evaluation  of  model  parameters 

There  are  three  commonly  used  methods  for  evaluating  model  parameters: 
the  moments  method  of  analysis,  the  frequency  response  analysis,  and  the 
s-plane  analysis.  The  advantages  and  disadvantages  of  the  moments  method 
have  been  discussed.   Comparison  of  the  frequency  response  analysis  and  the 
s-plane  analysis  has  been  given(14).  The  s-plane  analysis  is  applicable 
only  to  non-oscillatory  systems,  such  as  the  mixing  process  investigated  in 
this  work.   The  details  of  the  s-plane  analysis  have  been  described  by  the 
previous  worker  (i).  Model  parameters  were  evaluated  according  to  his 
procedure. 
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The  f-distribution  with  by-passing  across  the  P-mixing  unit  contains 
three  parameters,  T,,  3,  and  p  as  mentioned  in  Chapter  5.   The  dimensionless 
dead  time  could  be  determined  reasonably  accurately  from  the  recorded  res- 
ponse curves.   Dead  time  D  was  determined  by  taking  the  difference  between 
time  at  the  onset  on  the  response  of  the  output  Y(t)  and  that  of  the  input 
X(t).   The  dimensionless  dead  time,  T.,  was  then  calculated  by  dividing  D 
by  t,  the  mean  residence  time.   Tables  9  and  10  show  X ,  3,  and  p  for  the 
two  sets  of  data. 

Analysis  of  P  and  Q  values 

The  area  under  the  input  distribution  curve  should  theoretically  be 
equal  to  the  area  under  the  output  distribution  curve  for  every  experimental 
run.   In  the  course  of  this  investigation,  it  was  found  that  in  most  of  the 
runs,  the  areas  under  the  output  curves  were  larger  than  those  under  the 
input  curves.  The  first  set  of  data  indicates  that  the  average  value  of  the 
area  under  output  curves  was  10%  larger  than  the  average  value  of  the  areas 
under  the  input  curves,  and  15%  larger  for  the  second  set  of  data.   For  a 
closed  system  the  mean,  i/Q,  of  an  exit  age  distribution  or  the  difference  of 
the  means, ^Iiflvx,  between  the  output  distribution  and  input  distribution 
should  be  equal  to  one.   In  order  to  satisfy  these  requirements,  the  data 
adjustment  program  was  applied  to  the  raw  data  as  mentioned  previously.  As 
shown  in  Tables  5  and  6,  the  values  of  the  correction  factors,  P  and  Q,  are 
different  from  one  and  zero  respectively.  The  average  value  of  P  is  2.056 
for  the  first  set  of  data  and  2.847  for  the  second  set  of  data  (Tables  9  and 
10).  The  average  value  of  Q  is  0.294  for  the  first  set  of  data  and  0.366 
for  the  second  set  of  data.  The  average  values  of  P  and  Q  for  the  first  set 
are  closer  to  one  and  zero  respectively  than  those  for  the  second  set.  This 
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TABLE   9.      EXPERIMENTAL  VALUES  OF  THE   PARAMETERS,  T.,    0,  AND  p  OF  THE 

^-DISTRIBUTION  MODEL  WITH  BY-PASSING  ACROSS  THE  T-MIXING  UNIT 
FROM  RUN   1   TO  RUN  8. 


Run  Number 

3 

P 

**d 
Observed 

Error 
Square  x  10 

1 

.571 

18.40 

.678 

.125 

2 

.853 

5.521 

.521 

1.825 

3 

.797 

2.630 

.469 

.0995 

4 

.755 

3.265 

.469 

.619 

5 

.788 

3.263 

.521 

.0999 

6 

.736 

3.680 

.521 

.1401 

7 

.665 

6.691 

.521 

.3792 

8 

.878 

4.237 

.417 

.3442 

TABLE    10.      EXPERIMENTAL  VALUES  OF  THE  PARAMETERS,*"^,    0,  AND  p  OF  THE 

P-DISTRIBUTION  MODEL  WITH  BY-PASSING  ACROSS  THE P -MIXING   UNIT 
FROM  RUN   11   TO  RUN  18. 


Run  Number 

3 

P 

^d 
Observed 

Error 
Square  x  10 

11 

.553 

51.107 

.735 

416.850 

12 

.834 

16.268 

.630 

50.750 

13 

.658 

107.40 

.735 

20.395 

14 

.658 

3424.02 

.630 

4.881 

15 

.680 

7.17 

.630 

.192 

16 

.723 

126.11 

.682 

43.48 

17 

.674 

91.277 

.682 

43.48 

18 

.517 

1685.96 

.735 

23.03 
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may  bo  duo  partially  to  tho  smaller  difference  of  areas  between  input  and 
output  for  tho  first  set  (10%)  than  for  tho  second  set  (15%). 

Two  typical  experimental  runs,  one  from  the  first  set  and  the  other 
from  the  second  set  will  be  considered.   Fig.  25a  is  the  input  curve,  X(t), 
of  Run  1  (first  set)  which  is  redrawn  from  the  recorded  curve.   Fig.  25b  is 
tho  output  curve  of  Run  1  redrawn  from  the  recorded  curve.   The  area  under 
the  output  curve  is  U%   larger  than  that  under  the  input  curve.   In  order  to 
satisfy  the  two  requirements  as  mentioned,  the  correction  factors,  P  and  Q, 
were  computed.   The  points  located  on  the  figure  are  the  points  of  the  cor- 
rected response,  Y(0)  ,  for  P  equals  2.03  and  Q  equals  0.32.   Connection  of 
these  points  will  give  the  adjusted  output  distribution  curve.   The  adjusted 
output  curve  has  a  higher  peak  than  the  recorded  curve,  while  the  tail  of 
the  adjusted  curve  is  lower  than  the  recorded  curve.  The  values  of  the 
points  on  the  tail  of  the  adjusted  output  curve  are  very  small  and  even 
slightly  negative.  The  unexpected  negative  values  can  be  explained  by  Equa- 
tion (9-2)  in  which  if  the  correction  factor  Q  is  large  such  that  1-Q6  is 

negative,  then  Y(9)  can  be  negative. 

c 

Figure  26a  is  the  input  curve,  X(t),  of  Run  16  (the  second  set).   Fig. 
26b  is  the  output  curve,  Y(t),  of  Run  16.   The  area  under  the  output  curve 
is  12.8%  larger  than  that  under  the  input  curve.  The  correction  factor,  P, 
equals  2.9  and  Q  equals  0.36.   The  adjusted  curve  again  has  a  higher  peak 
and  lower  tail  than  those  of  the  recorded  curve.   The  magnitudes  of  the  nega- 
tive values  of  the  tail  of  Run  16  are  apparently  larger  than  those  of  Run  1. 
Those  negative  values  of  some  points  on  the  tail  may  contribute  an  appreciable 
negative  values  in  the  moments. 
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Analysis  of  the  moments 

The  numerical  values  of  the  mean,  ^/'qyx*  whicn  is  ecJual  theoretically 

tow  for  the  corrected  data,  is  very  close  to  one  as  expected.   The  numerical 

2  2 

values  of  ao         must  equal  <jA  of  the  exit  age  distribution  as  it  was  described 

v  YX  0 

2 
in  Chapter  5,  therefore,  a  positive  value  of  ^  aflYX  is  always  expected.  A 

2 
larger  value  of  a  0    represents  a  larger  spread  of  the  distribution  curve 

0  YX 

and  hence  represents  a  larger  extent  of  mixing.  The  first  set  of  data  shows 

2 
that  all  the  ^crflYX  are  positive,  while  the  second  set  of  data  is  negative  for 

most  of  the  experimental  runs.   The  negative  values  can  be  explained  by  the 

2 
fact  that  the  second  moment  about  the  mean  of  the  output,  0LV,  is  smaller  than 

6  Y 

2  2 

those  of  the  input  aAy»  The  reason  for  the  smaller  aAY  value  may  be  due  to 

the  low  degree  of  mixing,  and  also  due  to  the  negative  values  of  some  points 
which  appeared  on  the  tail  in  the  adjusted  output  curves. 

Let  us  analyze  in  more  detail  the  numerical  values  of  the  second  moments 
about  the  mean  for  each  experimental  run.   For  the  first  set  of  data,  Run  1 
through  Run  4  were  carried  out  separately  by  varying  the  position  of  the  out- 
let probe  over  the  cross-section  at  the  downstream,  while  the  inlet  probe  was 

fixed  at  the  center  of  the  cross-section  at  the  upstream.   Run  1  and  Run  2 

9 
have  values  of  a  at   of  0.0792  and  0.0669  which  are  slightly  different  from 
6YX 

2 
each  other.  The  values  of  ^<7AYX  for  Run  3  and  Run  4  are  very  close,  0.2052 

and  0.2075  respectively.   However,  comparing  Runs  1  and  2  with  Runs  3  and  4, 

2 
the  values  of  a  aflw  are  quite  different.   If  the  concentrations  across  the 

9YX 

2 
cross-section  are  very  uniform,  the  *cflw  values  for  each  run  should  be  close 

0  yx 

or  even  equal  to  each  other.   Run  5  through  Run  8  were  also  recorded  separately 

by  varying  the  position  of  the  inlet  probe  across  the  cross  section,  while  the 

2 
outlet  probe  was  fixed  at  the  center.  The  values  of  *aAw  ^or  Run  '  through 
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Run  8  arc  0.1418,  0.1608,  0.1587,  and  0.1233  respectively.   Those  values  are 

2 
all  slightly  different.   Taking  the  average  value  of  4 Onyx  for  Run  1  through 

Run  4  and  the  average  for  Run  5  through  Run  8,  the  first  has  the  average 

value  of  0.1397  and  the  second  0.1457.   After  taking  the  average  value,  the 

difference  in  the  second  moment  about  the  mean  for  the  average  of  Run  1 

through  Run  4  and  for  the  Run  5  through  Run  8  is  0.0062  which  is  so  small 


that  it  is  almost  negligible.   If  one  takes  the  overall  average  for  the  eight 

2 
9YX 


2 
runs  of  the  first  set  of  data,  the^aQVY  value  is  equal  to  0.1429. 


2 
The  result  may  represent  the  true  variance,  a  (7Aw»  at  tn*s  experimental 

condition  (excluding  the  difference  in  positions),  because  it  is  a  statistical 
average  of  the  variances  evaluated  under  different  detecting  positions.   How- 
ever, to  eliminate  the  stratification  effect  on  the  evaluation  of  the  true 
variance  by  this  procedure  may  require  excessive  time.  An  improved  tracer- 
monitoring  technique  specifically  used  for  this  non-ideal  sieve  tray  is 
necessary.   It  will  be  discussed  later. 

Among  the  second  set  of  experiments,  Runs  11  through  14  were  carried  out 
separately  by  varying  the  position  of  the  outlet  probe  over  the  cross  section 

at  the  downstream,  while  the  inlet  probe  was  fixed  at  the  center  of  the  cross 

2 
section  at  the  upstream.  The  values  of  An    for  these  runs  are  -0. 1843, 

9YX 

-0.0600,  -0.0348,  and  -0.0646  respectively.  The  average  value  is  -0.0859. 

Runs  15  through  18  were  recorded  separately  by  varying  the  position  of 
the  inlet  probe  across  the  cross  section  at  the  upstream,  while  the  outlet 

probe  was  fixed  at  the  center  of  the  downstream  cross  section.  The  values  of 

2 
Acrnw  for  these  ™ns  are  0.0906,  -0.0043,  -0.5280,  and  0.0161  respectively. 

2 
The  average  value  of  aq         is  -0. 1067. 

y  ia 

The  difference  between  the  average  variance  for  Run  11  through  14  and 
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that  for  the  Runs  15  through  18  is  -0.0208  which  can  no  longer  be  negligible 
at  the  low  degree  of  mixing  which  occurred  in  these  runs. 

Table  11  contains  the  values  of  the  average,  variance,  and  standard 

2  2 

deviation  for  a  aflW  and  model  parameters.  The  average  of  AoflW  for  the  first 
y  Y/v  0  YA 

set  is  0. 1429  and  for  the  second  set  is  -0.0962.  This  shows  that  the  mixing 
for  the  first  set  of  runs  was  closer  to  backmix  than  that  for  the  second  set 
of  runs.   Recall  that  the  first  set  of  runs  was  carried  out  at  a  higher  air 
flow  rate  than  the  second  set  of  runs.  The  number  of  holes  and  the  flow  rate 

of  air  should  contribute  significantly  to  the  extent  of  agitation  and  stirring 

2 
and  hence  change  the  degree  of  liquid  mixing.   The  negative  value  of  4q 

0  yx 

for  the  second  set  of  runs  may  be  due  to  the  stratification  which  was  pre- 
dominate at  a  low  degree  of  mixing.  As  mentioned  in  the  previous  chapter, 
the  stratification  downstream  was  clearly  observed  when  the  extent  of  agita- 
tion in  the  bubbling  zone  was  small. 

2 
The  variance  of  a<j    for  the  first  set  is  0.132  and  the  standard  devia- 

0  YX 

tion  is  0.3638;  while  for  the  second  set,  the  variance  is  5.122,  the  standard 
deviation  is  2.263,   It  is  obvious  that  the  second  set  of  data  has  a  much 
larger  variance  (or  standard  deviation)  than  that  of  the  first  set  of  data. 
This  may  mean  that  the  distribution  of  the  tracer  concentration  for  the 
second  set  is  not  as  uniform  as  the  first  set  of  runs.  The  nonuniformity 
of  tracer  concentration  is  expected,  because  a  thorough  agitation  or  stirring 
probably  can  not  be  obtained  by  using  only  one  row  of  holes  for  air  passage. 
From  the  above  argument,  one  may  conclude  that  the  stratified  elements 
observed  half  way  between  the  inlet  and  bubble  zone  are  not  completely  elimi- 
nated by  the  bubble  zone;  the  effect  of  stratification  on  the  degree  of  liquid 
mixing  still  exists.  The  effect  can  be  decreased  by  increasing  the  degree  of 
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agitation,  for  instance,  increasing  the  percentage  of  free  area. 

Analysis  of  the  model  parameters 

Two  different  distribution  curves  can  have  equal  values  of  the  variance. 
In  order  to  characterize  the  shape  of  the  distribution  curve  and  also  the 
physical  significance  of  the  distribution  curve  in  this  system,  an  exit  age 
distribution  function,  namely  the  P-distribution  model  with  by-passing  across 
the  P-mixing  unit  was  used  in  this  work  as  stated  previously. 

The  parameter,  3,  in  the  P-mixing  model  with  by-passing  across  the 
P-mixing  unit,  is  defined  as  the  ratio  of  the  volumetric  flow  rate  through 
the  P-mixing  unit  to  the  total  volumetric  flow  rate.  The  value  of  1-3  hence 
represents  the  ratio  of  the  volumetric  flow  rate  of  the  bypassing  stream  to 
the  total  volumetric  flow  rate.   In  other  words,  the  value  of  1-3  represents 
the  fraction  of  the  entering  fluid  being  bypassed.  The  average  value  of  3 
for  the  first  set  of  data  is  0.755,  i.e.  24.5%  of.  the  entering  fluid  being 
bypassed,  while  the  average  value  of  3  f°r  the  second  set  of  data  is  0.662, 
i.e.  33.8%  of  the  entering  fluid  being  bypassed.  The  variance  and  standard 
deviation  of  3  for  the  first  set  of  data  are  0.0175  and  0.1323,  while  those 
for  the  second  set  of  data  are  0.0213  and  0.1462  respectively,  (see  Table  11). 
It  appears  that  the  nonuniformity  of  the  tracer  concentration  distribution 
or  the  stratification  effect  has  only  slight  effect  on  the  bypassing  para- 
meter. 

The  parameter,  p,  of  the  P-mixing  unit  has  a  physical  meaning  related 
to  the  number  of  pools  connected  in  series.  The  parameter,  p,  can  be  any 
value  equal  to  or  greater  than  one,  but  not  necessarily  integer  in  this  model. 
The  physical  significance  constrains  that  p  must  be  greater  than  one.  The 
parameter  p  represents  complete  mixing  or  backmix  when  it  is  equal  to  one, 
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and  plug  flow  when  it  is  equal  to  infinity.   The  values  of  all  the  p's 
obtained  in  this  work  are  greater  than  one.   This  indicates  that  the  use  of 
the  P-distribution  model  with  bypassing  across  the  P-mixing  unit  can  pro- 
perly describe  the  system.   The  average  values  of  p  for  the  first  set  and 
the  second  set  are  5.96  and  688.66  respectively.   The  average  value  of  p  for 
the  first  set  of  data  is  much  smaller  than  that  of  the  second  set  of  data. 
This  again  means  that  the  mixing  conditions  for  the  first  set  of  runs  is 
closer  to  the  completely  mixed  state  while  that  for  the  second  set  of  runs  is 
closer  to  plug  flow.   Comparison  of  the  p  values  of  these  two  sets  of  runs 
enables  one  to  conclude  that  the  p  value  will  become  close  to  unity  as  the 
number  of  holes  increases. 

The  value  of  the  variance  of  p  for  the  first  set  is  0.779  and  that  for 
the  second  set  is  3»256.  The  standard  deviations  for  the  first  set  and  the 
second  set  are,  therefore,  0.8825  and  1.8045  respectively.  This  also  means 
that  the  stratification  affects  the  model  parameter,  p.   However,  the  effect 
of  stratification  on  the  degree  of  liquid  mixing  and  the  values  of  model 
parameter  p  may  be  decreased  by  increasing  the  number  of  holes. 

T,  represents  the  dimensionless  dead  time.  T.  was  defined  in  the  model 

d  q 

as  the  total  volume  of  the  plug  flow  unit  divided  by  the  volume  of  the  system. 
As  has  been  described,  the  dead  time  D  was  experimentally  determined  by  taking 
the  difference  between  times  at  the  incipient  responses  or  at  the  onsets  of 
the  responses  of  the  output  Y(t)  and  input  X(t).  The  dimensionless  dead 
time,  T,,  was  then  calculated  by  dividing  D  by  t,  the  mean  residence  time. 
Theoretically  for  the  plug  flow  model,  the  difference  between  time  at  the 
onset  on  the  responses  of  output  Y(t)  and  input  X(t)  will  be  the  time  required 
for  each  fluid  element  to  travel  from  the  inlet  to  the  outlet.  Therefore, 
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the  dead  time  is  equal  to  the  mean  residence  time,  and  the  dimensionless  dead 
time  will  hence  be  equal  to  unity.  As  a  backmix  model,  the  dimensionless 
dead  time  approaches  zero  because  there  is  no  concentration  gradient  exist- 
ing in  the  backmix  system.   The  reasonable  value  for  Td  hence  lies  between 
one  and  zero.   The  average  value  of  the  first  set  of  data  is  0.515,  and  0.682 

for  the  second  set  of  data  (see  Table  11).   These  values  also  show  that  "C. 

a 

of  the  first  set  is  smaller  than  that  of  the  second  set.   This  appears  to  be 

reasonable  because  the  first  set  of  runs  is  closer  to  backmix  than  the  second 

set  of  runs  as  mentioned  previously.  The  variance  of  t*.  for  the  first  set  is 

0.02175  and  0.00507  for  the  second  set.  The  corresponding  standard  deviations 

of  T.  for  the  first  set  is  0.1475  and  for  the  second  set  is  0.0712.  The  vari- 
d 

ances  (or  standard  deviations)  for  the  first  set  and  the  second  set  are  all 

small.   It  appears  that  the  effect  of  stratification  on  the  model  parameter 

T,  is  small, 
a 

Conclusions 

1.  The  pulse  testing  technique  using  sodium  chloride  solution  as  a 
tracer  material  has  been  shown  to  be  an  effective  experimental  technique  for 
gathering  information  on  liquid  mixing  in  a  flow  system.  A  tremendous  amount 
of  time  and  money  can  be  saved  by  using  the  conductivity  measurement  instead 
of  the  measurement  of  radioactivity,  which  was  used  in  the  previous  work. 
The  response  of  the  instrument  required  for  the  conductivity  measurement  can 
also  be  made  sufficiently  fast.  The  instrument  used  in  this  work  has  a  time 
constant  of  about  0.25  sec.  The  advantages  can  be  listed  as  follows: 

a.  Salt  tracer  is  much  cheaper  than  radioactive  isotopes. 

b.  The  source  of  salt  tracer  is  almost  unlimited. 

c.  Biological  effects  are  negligible. 
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d.   Waste  disposal  problems  can  be  avoided. 

2.  The  degree  of  liquid  mixing  and  the  parameters  in  the  mixing  model 
are  affected  by  a  stratification  effect,  when  the  experiment  is  carried  out 
using  only  a  pair  of  probes. 

3.  The  extent  of  the  stratification  effect  on  the  degree  of  liquid 
mixing  and  the  model  parameter,  p,  can  be  decreased  by  increasing  the  agita- 
tion and  stirring  of  the  system,  for  instance,  an  increase  in  number  of 
openings  or  in  the  size  of  free  area  for  vapor  phase  passage. 

4.  The  P-distribution  model  with  by-passing  across  the  P-mixing  unit 
adequately  describes  the  liquid  mixing  on  a  distillation  tray.   The  P-mixing 
model  with  by-passing  from  the  inlet  to  the  outlet  is  also  expected  to  be 
adequate  in  describing  liquid  mixing  on  a  full  size  sieve  plate  distillation 
column. 

5.  The  moments  method  of  analysis  is  a  convenient  way  for  evaluating 
the  degree  of  liquid  mixing. 

6.  The  s-plane  analysis  is  a  suitable  method  for  the  evaluation  of 
model  parameters. 

7.  The  tracer  monitoring  can  be  carried  out  by  inserting  several  probes 
over  the  same  cross  section  at  the  downstream  or  at  the  upstream.   The  con- 
centration-time curve  can  be  obtained  by  taking  the  average  concentration 
from  the  curves  (recorded  from  the  responses  of  the  probes)  at  the  corres- 
ponding time. 

Figures  27a, b,c,  and  d  show  the  hypothetical  response  curves  of  the  four 
probes  located  at  equal  distance  over  the  cross  section  of  the  liquid  stream. 
One  can  combine  these  four  response  curves  to  obtain  a  single  response  curve 
by  simply  taking  the  arithmetic  average  of  the  four  responses  at  various 
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times.   For  instance,  let  one  define  C(t)  as  a  single  response,  C,  (t),  C2(t), 

C.(t),  and  C. (t)  are  responses  for  Fig.  27a, b,c,  and  d  respectively.   One 
->         4 

has,  for  the  combined  response, 

C(0)     =      [C^O)   +  C2(0)   +  C3(0)   +  C4(0)]  x  i 

a     TO  +0+0+0]  xi 
*-  4 

=     0 

C(l)     =     [0  +  200  +  100  +  0]  x  £ 

=   7"5S 

C(10)     =     fl980  +  2100  +  20801  *  T 
=     2090      . 


HI 


Fig.  27  a 
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Fig.  27a,  b, c,d.    Hypotherical   response     curves 
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CHAPTER  X 

PROPOSAL  FOR  FUTURE  STUDY 

Further  study  of  liquid  mixing  on  distillation  trays 

A  correlation  chart  or  working  equations  for  the  reliable  evaluation 
of  the  degree  of  liquid  mixing,  which  is  an  important  factor  of  the  overall 
tray  efficiency  would  be  useful  to  designers.  This  chart,  or  equations, 
should  show  the  relationship  between  the  model  parameters  and  the  four 
variables,  i.e.  liquid  rate,  gas  rate,  weir  height,  and  percentage  of  free 
area.  It  would  be  convenient  to  prepare  a  correlation  chart  in  such  a  way 
that  the  relationship  between  the  parameters  and  the  dimensionless  groups, 
such  as  the  Reynolds  number  and  the  Froude  number,  etc.  are  shown.   In  order 
to  prepare  a  good  correlation  chart,  a  large  amount  of  additional  data  is 
required.  The  data  should  include  several  different  sources:  data  from  the 
previous  study  (1),  data  from  the  non-ideal  sieve  tray  used  in  this  work,  and 
data  from  a  multistage  distillation  column.   Comparison  of  the  results  of 
these  different  sources  of  data  can  provide  a  more  reliable  industrial  appli- 
cation of  the  chart* 

Additional  data  for  liquid  mixing  on  non-ideal  sieve  tray 

1.  Considerable  attention  should  be  given  to  setting  up  a  good  statis- 
tical design  for  experimental  runs  before  any  data  are  taken.  A  factorial 

design  (see  Appendls  3)  is  recommended  here.  This  factorial  design  should 

4 
include  four  levels  for  each  variable,  i.e.  4  =256  runs,  not  including  repli- 
cate runs.  We  might  start  with  a  three  level  design  for  each  variable  which 
requires  3  =81  runs  as  shown  in  the  shaded  area  of  Appendix  4. 

2.  The  ^-distribution  model  with  by-passing  across  the  P-mixing  unit 
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has  boon  proved  to  bo  adequate  in  describing  the  flow  pattern.   Another  form 
of  the  P-mixing  model  with  a  different  path  for  by-passing  —  the  P-distribu- 
tion  model  with  bypassing  from  inlet  to  outlet  can  also  be  used  to  fit  the 
experimental  data.   Comparison  of  the  fit  of  the  two  models  will  clarigy 
further  the  mechanism  of  bypassing. 

3.  Another  interesting  investigation  is  to  insert  several  baffle  plates 
into  the  bubble  zone  at  different  depths  without  touching  the  bottom  of  the 
tray.   Changes  in  the  residence  time  and  bypassing  may  thus  be  detected. 
Decreases  in  bypassing  and  mixing  in  the  direction  of  flow  are  known  to  be 
favorable  to  the  increase  of  Murphree  liquid  efficiency. 

4.  The  actual  cause  of  the  oscillation  phenomenon  in  the  distillation 
tray  still  remains  unknown.   It  should  be  investigated  further. 

5.  The  investigation  may  be  repeated  using  a  non-Newtonian  fluid* 

Additional  data  for  liquid  mixing  on  multi-stage  distillation  column 

A  three-stage  distillation  column  is  being  built  in  the  department  of 
chemical  engineering,  Kansas  State  University.   This  will  enable  one  to  study 
liquid  mixing  on  actual  distillation  trays  by  using  the  conductivity  method. 
The  study  may  include  sieve  trays,  bubble  cap  trays  and  valve  trays  etc. 

The  study  of  concentration  transients  in  response  to  flow  rate  changes 

The  study  of  concentration  transients  in  response  to  flow  rate  changes 
can  be  carried  out  in  a  tube  or  a  coil  using  experimental  techniques  described 
in  this  work.  A  tee-fitting  at  the  outlet  of  the  tube  separates  the  flow  into 
two  lines,  one  of  which  is  equipped  with  a  solenoid  valve  and  a  globe  valve, 
the  other  with  a  globe  valve.   Under  steady  state  condition,  the  concentration 
will  be  constant  at  a  fixed  point  if  the  effect  of  turbulent  fluctuation  can 
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be  neglected.   By  varying  the  settings  of  the  globe  valves  it  is  possible  to 
introduce  a  step  change  in  velocity  which  in  turn  leads  to  a  change  in  con- 
centration along  the  axial  direction  of  the  tube. 

If  an  electrolyte  is  used  as  a  tracer,  one  can  measure  the  transient 
concentration  response  by  using  several  well-designed  conductivity  probes 
located  at  different  positions  along  the  axial  direction  of  the  tube. 

Assuming  plug  flow,  no  axial  mixing,  perfect  radial  mixing,  constant 
physical  properties,  and  constant  overall  mass  transfer  coefficient  with 
length,  one  may  first  try  to  fit  the  following  equation  to  the  experimental 
data,  i.e.  (50) 


|f  +-[l+r(t)]!|    =    -  P0[l+br(t)]c 


where 


C  =  dimensionless  concentration 

t  =  dimensionless  time 

r(t)  =  function  of  flow  disturbance 

x  =  dimensionless  distance  along  the  tube 

b  =  rate  of  change  of  overall  mass  transfer  coefficient  with  fluid 

flow  rate . 

P  =  the  measure  of  effectiveness  of  mass  transfer,  which  equals 
o 

zero  if  there  is  no  mass  transfer  taking  place  from  the  wall 
(or  packing)  to  the  fluid. 

If  r(t)  is  a  step  change  in  flow,  r(t)  =  au(t),  the  above  equation  has 
an  exact  solution  as  follows  (50) 
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exp(-   P    [x  -  a(l    -   b)t")  )        ,        t     <     T^— 
o  L  J  1    +  a 

C      = 

,  rx(l   +ab)  -,  .  x 

exp(-   P    |  )  ,  t     >    - 

oLi+a-  1   -f  a 

where 

a  =  dimensionless  magnitude  of  step  change  in  flow  rate 

u(t)  =  unit  step  change 

The  effect  of  axial  dispersion  can  be  taken  into  account  by  modifying  the 
governing  equation  to 

2 
l£  +  E  (t)  1?  +   n  +r(t)]i£  =  -  P  [i  +br(t)]c 

where 

E  (t)  =  the  axial  eispersion  coefficient  and  is  a  function  of  time 
x 

in  this  case 
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Appendix  la.  Rotameter  I ,    Calibration   curve  (Sp.Gr.=  I  ) 
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APPENDIX  III 
DESIGN  OF  THE  CONDUCTIVITY  INSTRUMENT  SYSTEM 

The  conductivity  measurement  system  used  in  this  work  were  designed  by 
the  Industrial  Instruments  Inc.,  Cedar  Grove,  N.  J..   This  system  contains 
two  Type  RA4  Solu  Meters  and  two  probe  type  conductivity  cells. 

The  following  is  a  brief  description  of  the  instruments  and  the  prices: 

Type  RA4-WB  Solu  Meter,  range  0-2500  micromhos/cm 

specific  conductance  refer  to  25°  C  with 
meter  0-25  evenly  divided  and  0-10  rav  DC 
output  proportional  to  conductivity;  for 
115  V  60  cycle  AC  operation  at $  231.00  ea. 

Probe  type  conductivity  cells,  epoxy, 

approxime  tely  12"  long  and  5  mm  OD; 

platinum  electrodes  and  cell  constant 

of  0.50/ cm  at $  55.00  ea. 


120 


CO 

c 

3 

»_ 

__ 

o 

+- 

c 

CD 

£ 

i_ 

CD 

Q. 

X 

CD 

CD 

J/> 

+- 

CD 

o 

> 

4_ 

CD 

T) 

rO 

-J 

cr 

i_ 

o 

■i  ■■ 

V- 

o 

**- 

_l 

CD 

o 

CD 

CD 

u_ 

#* 

O 

■4- 

"O 

X- 

CD 

CD 

u> 

CT 

-o 

<1> 

O 

CD 

o 

-C 

C 
CD 

o 

i_ 

< ) 

•  — 

\_ 

(/> 

•  • 

5 

CD 

Q. 

o 

UJ 

X 

Li. 

;> 

£ 

0. 

o 

2 

o 
u 

Gq;    ;o    uSisap   idojjsiws    pasodoJd  'trXipuaddv 


121 


ACKNOWLEDGEMENT 

The  author  gratefully  acknowledge  the  able  direction  and  inspiring 
guidance  of  his  thesis  adviser,  Dr.  Liang-tseng  Fan. 

He  thanks  Dr.  John  L.  Johnson  for  his  suggestions  and  help  in  preparing 
the  conductivity  instrument. 

A  deep  debt  of  gratitude  is  also  due  to  the  members  of  the  supervisory 
commettee,  Dr.  William  H.  Honstead,  Head,  Department  of  Chemical  Engineering, 
Kansas  State  University,  Dr.  John  C»  Matthews,  and  Dr.  Lyle  J.  Dixon,  for 
their  encouragement  and  help. 

Grateful  acknowledgement  is  also  due  to  Dr.  Richard  G.  Akins  for  his 
suggestion  with  part  of  the  equipment  set  up. 

Financial  assistance  for  this  thesis  was  provided  by  the  National 
Science  Foundation  under  NSF  grant  GK67.  The  author  is  grateful  for  this 
generous  support. 


122 


LITERATURE  REVIEW 


1.  Johnson,  J.  L. ,  Ph.D  Thesis,  Department  of  Chemical  Engineering,  Kansas 
State  University,  Manhattan  (1966). 

2.  Danckwerts,  P.  V.,  Chem.  Eng.  Sci.,  1_,  1  (1953). 

3.  Noar,  P.  and  R.  Shinnar,  I&EC  Fundamentals,  2,  4,  (1963). 

4.  Bell,  R.  L.  and  A.  L.  Balb,  On  the  interpretation  of  measured  residence 
time  distributions  in  chemical  process  equipment.   A.I.Ch.E.  Reprint 
22E  Fifty-eight  annual  meeting,  December  5-9  (1965). 

5.  Danckwerts,  P.  V.,  Chem.  Eng.  Sci.  2  (1958). 

6.  Spalding,  D.  B. ,  Chem.  Eng.  Sci.  9,  74  (1958). 

7.  Zweitering,  T.  N. ,  Chem.  Eng.  Sci  H_,  1  (1959). 

8.  Levenspiel,  0.  "Chemical  Reaction  Engineering,  John  Wiley  &  Sons  Inc., 
New  York  (1962). 

9.  Eguchi,  W. ,  Proc.  25th  Anniv.  Congress,  Soc.  Chem.  Eng.  (Japan),  Nov.. 
(1961). 

10.  Sater,  V.  E.  and  Octave  Levenspiel,  I&EC  Fundamentals  5_,  1  (1960). 

11.  Leonard,  E.  F. ,  Ph.D  dissertation,  Department  of  Chemical  Engineering, 
University  of  Pennsylvania,  Philadelphia  (1960). 

12.  Hay,  J.  R. ,  M.S.  Thesis,  Department  of  Chemical  Engineering,  Vandervilt 
University,  Nashville,  Tenn  (1964). 

13.  Hay,  J.  R. ,  A.I.Ch.E.  Annual  meeting,  Philadelphia  Dec  5-9  (1965). 

14.  Clements,  W.  C.  Jr.,  Ph.D  dissertation,  Vandervilt  University,  Nashville, 
Tenn.  (1963). 

15.  Otto,  R.  E.  and  Jr.  Stout,  Chemical  Engineering  Progress  Symposium 
Series,  Process  Dynamics  and  Control,  36,  57  (1961). 

16.  Moser,  J.  H.  and  C.  R.  Cupit,  Fifty-eight  National  Meeting,  Reprint  5D 
Dallas,  Texas,  Feb.  6-9  (1966). 

17.  Bishoff,  K.  B.  and  E.  A.  Mccraken,  Principles  of  the  interpretation  of 
tracer  tests  of  flow  systems  part  1  &  2. 

Manuscript  from  I&EC  Research  Results,  Service  treat  as  a  private  com- 
munication (1966). 

18.  Adler,  R.  J.,  "Finite  Stage  Modeling".  Presented  at  the  Fourth  National 
Chemical  and  Petroleum  Symposium  of  the  Instr.  Society  of  America,  May 
(1961). 


123 


19.  Kirchnbum,  E. ,  "Distillation  and  Roctif icat ion",  Trans,  by  M.  Wulfinghoff, 
Brooklyn  Pub.  Co.  (1948). 

20.  Nord,  M.,  Trans.  Am.  Inst.  Chem.  Engrs.,  42_,  863  (1946). 

21.  Gatrcaux,  M.  F. ,  and  H.  E.  O'Connoll,  Chem.  Eng.  Prog.,  5J_,  232  (1961). 

22.  Micklcy,  H.  S.,  L.  A.  Gould,  and  L.  M.  Schwartz,  Canadian  Journal  of 
Chemical  Engineering.,  137,  14  (1961). 

23.  Oliver,  E.  D.,  and  J.  Marangozis,  Can.  J.  of  Chem.  Eng.,  36,  161  (1958). 

24.  Warzel,  L.  A.,  Ph.D  Dissertation,  University  of  Michigan  (1955). 

25.  Johnson,  A.  I.,  and  J.  Marangozis,  Can.  J.  of  Chem.  Eng.,  3_6,  161  (1958). 

26.  Anderson,  J.  E. ,  Sc.D.  Dissertation,  M.I.T.  (1954). 

27.  Gerster,  J.,  et.al.,  "Tray  Efficiency  in  Distillation  Column",  Final 
Report  from  University  of  Delaware,  A.I.Ch.E.,  New  York  (1958). 

28.  Welch,  N.  E.,  Ph.D  Dissertation,  Texas  A  &  M  University,  College  Station 
(1958). 

29.  Iyer,  S.  R.  and  P.  S.  Murti,  Indian  J.  Technol.  3,  75  (1964). 

30.  Foss,  A.  S.  et.al.  A.I.Ch.E.  Journal  4,  231  (1958). 

31.  Mutzenberg,  A.,  Chem.  Eng.  Techn.  34,  542  (1962). 

32.  Olson,  R.  E.  and  W.  J.  Heidiger,  University  of  Washington,  Seattle  (1962). 

33.  Strand,  C.  P.,  Chem.  Eng.  Prog.  59,  58  (1963). 

34.  A.I.Ch.  E.  "Tray  Efficiencies  in  Distillation  Columns,"  Final  Reports  of 
Research  Committee  (1958,  1959,  1960). 

35.  Cha,  L.  C. ,  Ph.D  Dissertation,  Department  of  Chemical  Engineering,  Kansas 
State  University,  Manhattan  (1965). 

36.  Perry,  J.  H. ,  R.  H.  Perry,  C.  H.  Chilton,  and  S.  H.  Kirkpatrick, 
"Perry's  Chemical  Engineers'  Handbook."  Fourth  Edition,  McGRAW-HILL  BOOK 
COMPANY  INC.,  New  York  (1963). 

37.  A.I.Ch.E.  Research  Committee,  "Bubble  Tray  Design  Manual"  (1958). 

38.  Brown,  G.  G. ,  "Unit  Operations"  Second  Printing,  John  Wiley  &  Sons  Inc., 
New  York  (1951). 

39.  Yakubovich,  V.  I.  Ulanov,  and  A.  V.  Machinskii,  Zavodsk.  Lab.,  29  (9), 
1141-3  (1963). 


124 


40.  Westerterp,  K.  R.  and  P.  Landsman,  Chem.  Eng.  Sci.  _17_,  4  (1962). 

41.  Woong  Ki  Kang  (At.  Energy  Res.  Inst.,  Seoul)  Ta  Han  Hsueh  Hui-Chih, 
6(2),  130-2  (1962). 

42.  Stenhagen,  E.  and  0.  Mellander,  Acta  Chem.  Scan.  1_0,  1317-26  (1956). 

43.  Prausnitz,  J.  M. ,  Ph.D  Dissertation,  Department  of  Chemical  Engineering, 
Princeton  University,  New  Jersey  (1955). 

44.  Manning,  F.  S.,  Ph.D  Dissertation,  Department  of  Chemical  Engineering, 
Princeton  University  (1959). 

45.  Schlichting,  H. ,  "Boundary  Later  Theory"  Fourth  Edition,  McGRAW-HILL 
BOOK  COMPANY,  INC.,  New  York  (1960). 

46.  24th  Meeting  of  the  European  Federation  of  Chemical  Engineering  Proceed- 
ings of  the  International  Symposium  on  Distillation,  Brighton,  England 
4/5/6  May  (1960). 

47.  Teller,  A.  J.  and  R.  E.  Rood,  A.I.Ch.E.  Journal,  July  (1962). 

48.  Cheng,  S.  I.,  and  A.  J.  Teller,  A.I.Ch.E.  Journal,  7,  282  (1961). 

49.  Rooze,  J.  M. ,  Bulletin  SRC  20-C-62-6,  Systems  Research  Center,  Case 
Institute  of  Technology  (1962). 

50.  Kamman,  D.  T.  and  L.  B.  Koppel,  I&EC  5,  2  (1966). 

51.  MCALLISTER,  R.  A.  and  C.  A.  Plamk,  A.I.Ch.E.  Journal  4,  3  (1958). 


AN  EXPERIMENTAL  INVESTIGATION  OF  LIQUID  MIXING 
ON  A  NON IDEAL  SIEVE  TRAY 


by 


TSUNG-WEN  LEE 
B.  S.,  National  Taiwan  University,  Taiwan,  China,  1962 


AN  ABSTRACT  OF  A  MASTER'S  THESIS 

submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree 


MASTER  OF  SCIENCE 


Department  of  Chemical  Engineering 


KANSAS  STATE  UNIVERSITY 
Manhattan,  Kansas 


1967 


A  12  x  48"  rectangular  sieve  tray  column  with  rectangular  downcomers 
was  used  to  study  liquid  mixing  effects  on  a  non-ideal  plate.   Using  the 
air-water  system,  liquid  mixing  tests  were  made  by  using  the  pulse  testing 
technique.  Sodium  chloride  was  used  as  a  tracer.  Data  were  taken  for  a 
tray  with  one  row  of  holes  (15  holes)  and  three  rows  of  holes  (43  holes). 

Various  flow  patterns  were  visually  observed  and  were  also  photographi- 
cally recorded  and  analyzed.  By  using  conductivity  measurements,  the  effect 
of  stratification  on  the  degree  of  liquid  mixing  was  investigated.   Using 
the  technique,  a  tracer  monitoring  technique  was  developed.  Experimental 
data  were  analyzed  by  using  the  moments  method  in  which  the  values  of  the 

difference  of  the  second  moments  about  the  mean  between  the  input  and  output 

2 
pulses,  afl,  were  calculated.  The  V -distribution  model  with  bypassing  and 

9 

a  modified  model  with  a  different  path  for  bypassing  were  employed  to  fit 
the  experimental  data. 

Future  work  has  been  proposed,  which  includes  a  statistical  design  of 
experimental  work. 


